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Introduction 



Basis of the general theory of differential equations are the theory of solutions and the 
theory of integrals. The functional-analytical research of integrals is most deeply developed 
for ordinary differential systems and linear systems of partial differential equations. 

The initial problem of the integration in quadratures for differential equations has led to 
necessity to develop methods of analytical and qualitative researches both for solutions and 
integrals of differential systems. So, J. Liouville [4, 5] considered the problem of the integration 
in quadratures for the Riccati equation. His investigation gave the classical problems about 
a form of solutions and about development of methods for finding solutions of given forms. 

For the first time the problem of building a general integral from first integrals was con- 
sidered by J. Jacobi in [6-8]. He also introduced the notion of a last multiplier (also known 
in publications as Jacobi's last multiplier) and he used this notion of a last multiplier to solve 
the problem of finding a general integral. 

The profound researches, which are the base of the theory of integrals, are due to F.G. 
Minding [9], A.V. Letnikov [10], and A.N. Korkine [11]. 

One of such problems is the Darboux problem [12] about building of a general integral 
by known partial integrals and about the form of a general integral for the ordinary differ- 
ential equation of the first order. The Darboux problem for ordinary differential systems, 
total differential systems, and systems of partial differential equations was considered in the 
monographs [1, 2, 13, 14] and in the articles [15 - 46]. 

The method of last multiplier was developed by S. Lie. He has given the new interpretation 
of this method and has created the theory of infinitesimal transformations in [47, 48]. In his 
papers were selected the fundamental approaches of the integration. These approaches are the 
base of the theory of closed differential systems. Methods of the integration were considered 
in [49] with regard to the uniform positions of the group analysis. This approach is gave the 
theoretical base for the classification [50, 51] of these methods. 

At the beginning of the twentieth century the interest to global researches for differential 
systems has decreased. But in the middle of the twentieth century the interest in it has 
appeared again. Let us note only the monographies: N.M. Gjunter [52], E. Cartan [53, 54], 
N.G. Chebotarev [55], L. Eisenhart [56], P.K. Rashevskii [57], H. Cartan [58], A.S. Galiullin 
[59], L.V. Ovsiannikov [60], N.P. Erugin [61], E.A. Barbashin [62], A.M. Samoilenko [63], 
P. Olver [64], A. Goriely [65]. It happened for the reason that these results find applications 
in the mathematical physics [66 - 71]. 

The intensive development of the qualitative theory of differential equations (the base of 
this theory was founded by H. Poincare [72]) also has led to solving some problems for the 
theory of integrals. So, the analytic structure of integrals and of integrating multipliers in a 
neighbourhood of a center was obtained by A.M. Liapunov [73] and N.A. Sakharnikov [74] 
respectively. Investigation of limit cycles for differential systems on the base of partial integrals 
and integrating multipliers was considered by M.V. DoIoa' [75 - 77]. Behaviour of trajectories 
for ordinary autonomous differential systems of the second order having the integral curves of 
the concrete forms and with the special qualitative properties was researched in [16; 78 - 90]. 

The subject of our investigation is a system of total differential equations 



where t G W^, x G M", m ^ n, dt = colon (dt^ , • • • , dt^^, dx = colon (dx^, . . . , dx^), X{t, x) = 

= (t, x) \\ is an n x m matrix with entries X^j : H R, i = 1, . . . ,n, j = 1, . . . ,m, H 

is a domain of the extended space M™+"; a linear homogeneous system of partial differential 
equations 



dx = X{t, x) dt, 



(TD) 



2.{x)y = 0, j = l,...,m, 



id) 



where y G M, x G M", m ^ n, the linear differential operators of first order 
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S^j{x) = ^ Uji{x) for all X G G, j = 1, . . . ,m, (0.1) 

1=1 

with coordinates Uj^: G W, j = 1, . . . ,m, i = 1, . . . ,n, a domain G C M"; and a Pfaff 
system of equations 

LUjix) = 0, j = l,...,m, (Pf) 
where x G M", m ^ n, the linear differential forms 

n 

u;j{x) = ''^Wjj^{x) dx^ for all x € G, j = l,...,m, (0.2) 

i=l 

with coefficients lo^j : G ^ M, j' = 1, . . . , m, i = 1, . . . , ra, a domain G C R". 
We recall that by domain we mean open arcwise connected set. 

Let the linear differential operators (0.1) of the system of partial differential equations (d) 
and the linear differential forms (0.2) of the Pfaff system of equations (Pf) be not linearly 
bound on the domain G [91, pp. 105 - 115]. Note that the operators , j = l,...,m 
(the 1-forms ujj , j = 1, . . . ,m) are called linearly bound on the domain G if these operators 
(1-forms) are linearly dependent in any point of the domain G [60, pp. 113 - 114]. 

The system (TD) is induced the m linear differential operators of first order 

n 

Xj{t,x) = dt,+^X^j{t,x)d^, for all (t, x) G n, j = l,...,m. (0.3) 

i=l 

We'll say that these operators are operators of differentiation by virtue of system (TD). 

Under the condition m = 1, we have the system (TD) is an ordinary differential system 
of n order. 

The system (TD) is said to be completely solvable on a domain 11' c 11 if for any point 
(tg ,Xq) G n' there exists a unique solution of the Cauchy problem with initial data (tg ,Xq) 
[1, p. 17]. If n' = n, then we say tliat the system (TD) is completely solvable. 

Suppose X G G^(n), i.e., the functions X- , i = 1, . . . , n, j = 1, . . . , m, are continuously 
differentiable on the domain H. Then the system (TD) is completely solvable if the Frobenius 
theorem [1, pp. 17 - 25; 58, pp. 290 - 297; 92, pp. 309 - 311; 93, p. 21] is true. 

Theorem 0.1(ti3e Frobenius theorem). The system (TD) with X G G^(n) is completely 
solvable if and only if the Frobenius conditions hold 

n n 

d^Xi^it, x) + Yl ^)^x,^ic(*' ^) = x) + Y^ %(i, x)d^X,^{t, x) 

(0.4) 

for all (t, .t) G n, i = 1, . . . ,n, j = 1, . . . ,m, = 1, . . . ,m. 

Using the operators (0.3), we obtain the Frobenius conditions (0.4) are represented via 
Poisson brackets as the system of the operator identities 

[Xj{t,x),X^{t,x)] =0 for all (t, x) G n, j = 1, . . . ,m, ( = 1, . . . ,m, (0.5) 

where O is the null operator. 

The system (TD) is the Pfaff system of equations 

Vi{t,x) = 0, i = l,...,n, (0.6) 

with the linear differential forms 

m 

?7j(t,x) = — ^Xj^(t,x) c/tj for all {t,x) eH, i = l,...,n. (0.7) 
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The Probenius conditions (0.4) for completely solvability of system (TD) (the Pfaff system 
of equations (0.6)) are represented via differential 1-forms (0.7) as the system of exterior 
differential identities [58, pp. 290 - 297] 

dr]^{t,x) A (^J\^r]^{t,x)j =0 for all (t, x) G n, i = l,...,n. (0.8) 

Consider an autonomous total differential system 

dx = X{x) dt, (ATD) 

where t G M™, x G M", m ^ n, dt = colon (dti, ... ,dt^), dx = colon (dx^, ... ,dxj^), the 
entries of the nx m matrix X{x) = \\X^j{x)\\ are X^j : G — >■ M, i = 1, . . . , n, j = 1, . . . , m, G 
is a domain of the phase space R". We say that the linear differential operators of first order 

n 

Xj{t,x) = dt,+J2^iji^)9x, for all (t,x) gR'" X G, j = l,...,m, (0.9) 

i=l 

induced by this system are nonautonomous operators of differentiation by virtue of system 
(ATD). The linear differential operators of first order 

n 

^^jix) = ^Xij{x)d^, forallxGG, j = l,...,m, (0.10) 

i=l 

are called autonomous operators of differentiation by virtue of system (ATD). 

Let XgC^(G). Then the Probenius conditions (0.5) for completely solvability of system 
(ATD) are the identities [1, pp. 112 - 113] 

[2lj(x),2l^(x)] =D forallxGG, j = l,...,m, C = l,...,m. (0.11) 

The system (d) with 2,j G C^(G), j = 1, . . . ,m (i.e., the coordinates Uj^ , j = 1, . . . ,m, 
i = l,...,n, are continuously differentiable on the domain G) is said to be complete [92, 
p. 521; 53, p. 117] if Poisson brackets of any two differential operators (0.1) are the linear 
combination of operators (0.1) 

m 

[£,(x),£^(x)] =^A^.^,(x)£,(x) forall xGG, j = l,...,m, C = l,...,m, (0.12) 
i/=i 

where the coefficients Aj^^^ G C^(G), j = 1, . . . ,m, ( = 1, . . . ,m, v = 1, . . . ,m. 
If Poisson brackets of operators (0.1) are symmetric, i.e., 

[£^ix),2^ix)] = [£^ix),2^{x)] for all xGG, j = 1, . . . ,m, C = I, ■ ■ ■ ,rn, (0.13) 

then we say that the linear homogeneous system of partial differential equations (d) with 
£j G G^(G), j = l,...,m, is jacobian [52, p. 62; 92, p. 523]. 

The symmetry (0.13) of the Poisson brackets of operators (0.1) is equivalent to 

[£j{x),£^{x)] =D forallxGG, j = l,...,m, C = l,---,m- (0.14) 

The identities (0.14) are the identities (0.12) with Aj^^{x) = for all x E G, j = 1, . . . ,m, 
C = 1, . . . , m, u = 1, . . . ,m. Therefore the jacobian system {d) is complete [52, p. 62]. 
A differential system 

d^,y = mj{x)y, j = l,...,m, (NS) 

where y G R, x G R", m ^ n, the linear differential operators of first order 
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Mj{x) = Uj^{x) d^^ for all X G G, j = l,...,m, (0.15) 

s=m+l 

is called a normal lineai homogeneous system of partial differential equations [52, p. 64]. 

A complete normal partial system is jacobian [52, p. 65; 1, pp. 38 - 40]. 

Problems of the theory of integrals for total differential systems were considered in process 
of necessity at the decision of adjacent problems. There are first of all the problems of orbits 
topology for completely solvable autonomous total differential systems were studied [93, 94]. 
Directly problems of the theory of integrals for multidimensional differential systems (TD), 
(d), and (Pf) are considered in [1, 2]. 

1. First integrals of total differential system 
1.1. First integral 

Suppose the system (TD) has the matrix X G G(n), i.e., the entries X- , i = 1, . . . ,n, 
j = 1, ... ,171, of the matrix X are continuous functions on the domain 11. 

Definition 1.1. A scalar function F G C^(n') is said to be a first integral on a domain 
n' C n of system (TD) with X G G(n) if the differential of the function F by virtue of 
system (TD) vanishes on the domain 11': 

dF{t,x)^^^^^ = for all {t,x) G W. (1.1) 
The differential of the function F by virtue of system (TD) is 

m n 

^^^^'^^[(TD) = E d,F{t,x)dt^ + J2d./it,x)dx,^^^^^ = 

7 = 1 1=1 




where Xj , j = 1, . . . ,m, are the operators (0.3). From (1.1) it follows that 

XjF{t,x) = for all (i,x) G n', j = l,...,m. (1.2) 

The connection between the identity (1.1) and the system of identities (1.2) is the jus- 
tification of that the linear differential operators (0.3) have been named the operators of 
differentiation by virtue of system (TD). 

Example 1.1. The total differential system 

dx^ = {x-^t^^ + ij^Xg) dti + t^X2 dt2, 

dx2 = {-I - Xjtj"^ + xl f^"^ + X2) dti + ( - 1 - x-^f^^ + xlf^"^ + x| - X2X3) dt2, (1.3) 

dx^ = X2{x.^ dt^+ (.X2 + Xg) dt2) 

induces the linear differential operators of first order 

%^{t, x) = + {x^t^^ + 1^X2) d^^ + {-l- x^t^^ + xl t^^ + xl) d^^ + X2X3 d^^ (1.4) 
and 

X2{t,x) = dt^ + 1^X28^^ + {-l-x^f^^ + xlt'{'^ + xl-X2X^)d^^ + X2{x2 + x^)d^^ (1.5) 
on the set D = {{t, x) : / 0} C M^. 
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The operations of operators (1.4) and (1.5) on the scalar function 

F: {t,x) ^ [l-xlt];^ -xl-xl)e^p{-2x^ti^) for all {t,x) e D (1.6) 

are identically equal to zero on the set D : XiF(t,x) = X2F{t,x) = for all {t,x) G D. 

By definition 1.1, the function (1.6) is a first integral on any domain 11 C D of the total 
differential system (1.3). 

The direct corollary of Definition 1.1 is the criterion of the existence of a first integral for 
a completely solvable total differential system. 

Theorem 1.1. A scalar function F G C^(n') is a first integral on the domain 11' C 11 
of the completely solvable on the domain H' system (TD) with X G C{Il) if and only if 
the function 11' — > M is constant along any solution x: t — > x{t) for all t G T' of 
system (TD), where the domain T' C M™' is such that {t,x{t)) G 11' for all t G T', i.e., 
F{t,xit)) = C for all t G T', C = const. 

If the system (TD) is not completely solvable, then the system (TD) can have first integrals 
even in the case when the system (TD) does not have solutions. 

Example 1.2. The autonomous linear system of total differential equations 

dxi = x-^ d,t-^ + d.t2 , dx2 = {1 + x^ + 2x2) dt^ + {x^ + 8x2) dt2 (1.7) 
has the linear differential operators of first order 

Xi(t, x) = + 2li(x) for all {t, x) G M^, X2{t, x) = d^^ + 2l2(x) for all {t, x) G M^, 
where 

2li(x) = x^ d^^ + (1 + Xi + 2x2) d^^ for all x G R^, 

2t2(x) = 3xi d^^ + (xi + 8x2) d^^ for all x G 

The Poisson bracket for the autonomous operators of differentiation by virtue of the total 
differential system (1.7) 

[2li (x), 2I2 (x)] = (8 - Xi) d^^ for all x G 

is not the null operator on any two-dimensional domain from the phase plane R^. By the 
Frobenius theorem (Theorem 0.1), the system (1.7) is not completely solvable. 

The operations of the nonautonomous operators of differentiation by virtue of system (1.7) 
on the holomorphic scalar function 

F: (i,x) Xiexp( - (ii + 8*2)) for all {t,x) G (1.8) 

are identically equal to zero on the space M^: XiF(t,x)= X2F{t,x)= for all {t,x) G M^. 
By definition 1.1, the function (1.8) is a first integral on the space of system (1.7). 
Let us prove that the system (1.7) has no solutions. 

Since x^: t ^ Ccxp(t^ + 8^2) for all t G M^, we see that the first equation of system 
(1.7) is an identity on the plane and the second equation of system (1.7) is 

dx2 = P{t,X2) dti + Q{t,X2) dt2 , (1.9) 

where 

P: {t, X2) ^ 1 + 2x2 + C'exp(ii + 8*3) for all {t, X2) G M^, 

Q: (t, X2) ^ 8x2 + C exp(ti + 8*2) for all {t, X2) G M^, 
C is a constant from the field M. 
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If the equation (1.9) possess the solution - t ^ X2i^) fo'^ i G T, where T is some 
domain of the plane M?, then the following conditions hold 

x^it) = P{t, x^it)) for all t G T, d^^ x^it) = Q{t, x^it)) for all t G T. 

But such the function T ^ M is not exist. It follows that the mixed derivatives of 
the second order 

\t, ^2(0 = \ Pit, X2it)) = (i) + 5Cexp(ti + 3*2) for all t G T, 

^t^t, ^2(*) = \ Qit, ^2(0) = 3 + 6x2 (t) + 4Cexp(ti + 3*2) for all t G T 

are not coincide neither at any C from the field M nor on any domain of the plane M^. 

Thus the not completely solvable total differential system (1.7) has the first integral (1.8), 
but the system (1.7) doesn't have solutions. 

1.2. Basis of first integrals 

Consider the set of scalar functions 

F, ■ {t, x) ^ Fjt, x) for all (t, x) G B', C C\n), s = l,...,k, B' C H C M™+", (1.10) 
and form the vector function 

F: {t,x) ^ {Fi{t,x),...,F^{t,x)) for all (t, x) G B' (1.11) 
with range EF cR'^. 

Theorem 1.2. Suppose the functions (1.10) are first integrals on the domain B' C B of 
system (TD) with X G C(B). Then the function 

^■.it,x)^<P{F^{t,x),...,F^{t,x)) foraU{t,x)eU', (1.12) 

where arbitrary scalar function $ G C^(EF), is a first integral on the domain B' of the total 
differential system, (TD). 

Proof. From Definition 1.1 for the first integrals (1.10), we have 

Xj F^{t, x) = for all (t, x) G B', j = 1, . . . , m, s = l,...,k. 

Then, for any scalar function $GC^(Ei^) on the range EiF of the vector function (1.11), 
we obtain 

k 

^{F{t,x)) = J2 dp, *(^)|^^^(,^,) % Psit,^) = for all {t,x) G B', j = l,...,m. 

s=l 

By Definition 1.1, the function (1.12) is a first integral on the domain B' of system (TD).B 

This theorem expresses functional ambiguity of a first integral for a system of total differ- 
ential equations: if a scalar function Fj^gC^(B') is a first integral on the domain B'cB of 
system (TD) with X G C(B), then the function {t,x) $(Fi(t,x)) for all {t,x) G B', 
where arbitrary scalar function $gC\EFi), is a first integral on the domain B' of system (TD). 

Thus the priority of first integrals functionally independent on some domain is installed. 
Therefore for a system of total differential equations we have the problems about the existence 
and the number of functionally independent first integrals. 

Definition 1.2. A set of the functionally independent on the domain B' C B first 
integrals (1-10) of system (TD) with X G C(B) is called a basis of first integrals on the 
domain B' of system (TD) if for any first integral on the domain B' of system (TD), 
we have ^(t,x) = ^(Fi{t,x), . . . ,Fi^{t,x)) for all (i,x)GB', where $ is some function of 
class C^(EF), EF is the range of the vector function (1.11). The number k is said to be the 
dimension of basis of first integrals on the domain B' of system (TD). 

A basis of first integrals we'll name also as an integral basis. 
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Definition 1.3. We'll say that two systems of total differential equations are integrally 
equivalent on some domain if on this domain each first integral of the first system is a first 
integral of the second system, and on the contrary each first integral of the second system is a 
first integral of the fi,rst system. 

The integrally equivalent on the domain 11' total differential systems have the same 
integral basis on this domain. 

1.3. Dimension of basis of first integrals for completely solvable systems 

Suppose the system (TD) has the matrix X G C°°(n). i.e., the entries X^ - , i = 1, . . . ,n, 
j = 1, . . . ,m, of the matrix X are liolomorpiiic functions on the domain 11. T is a domain 
from the space M"* such that for any solution x: t ^ x{t) for all t E T of the completely 
solvable system (TD), it follows that {t,x{t)) G IT for all t G T. Then, by the Cauchy 
theorem (see, for example, [93, p. 26]), for any point tQ G T there exists a neighbourhood 
Uq such that the solution x: t ^ x(t) for all t £\Jq of the completely solvable system (TD) 
with X G C°°(n) is a holomorphic function. Moreover, solutions of the completely solvable 
total differential system (TD) with X G C°°(n) depends holomorphically on parameters and 
initial data [93, pp. 39 - 41]. 

By x: t — > x{t; {tg ,Xq)) for all t G T we denote the solution x: t ^ x{t) for all i G T 
satisfying the initial condition x{tQ) = Xq. 

Lemma 1.1. Let the vector function x: t ^ x{t; {tg ,Xq)) for all t,tQ G T' be a solution 
on some simply connected domain T' C T of the completely solvable system (TD) with 
X G C°°(n). Then for any t^, G T' there exists the solution x: t x{t; (t^ , x^)) for all 
t £ T' of system (TD), where x^ = x{t^: {Iq ,.Xq)), such that x{tQ; (t^ ■,x{t^; {t-Q ,Xq)))) = Xq . 

Proof. Let x, x be the solutions of the corresponding Cauchy's problems of system (TD): 



The values x{t^) = x{t^). Then, by the Cauchy theorem, we have x{t) = x{t) for all t G T', 
where a simply connected domain T' such that the points ^^^^ belongs to T'. 
Thus x{t(i) = x{tQ) = Xq, i.e., x{tQ; {t^ ,x{t^; {t^ ,Xq)))) = Xq. M 

Lemma 1.2. Suppose the completely solvable system (TD) with X £0^(11) in a neighbour- 
hood of the point (tQ,.TQ)Gn satisfi,es the conditions of the Cauchy theorem. Then this system 
has n functionally independent on some neighbourhood of the point (tg , Xq) first integrals. 

Proof. Let x: t ^ x{t; (tg ,Xq)) for all t E T' be a solution of system (TD) on a simply 
connected domain T' 3 tg . The function 



x: t ^ x(t; (tp , Xq)) for all t eT' 



and 



X : t -> x{t; {t^ , x{t^; (tg , Xq)))) for all t G T'. 



F: {t,x) 



x{tQ\ (t, x)) for all (t,x) G U, 



00 ' 



where Uqq is some neighbourhood of the point (igjXg), under 



oo 




for all t G U, 



' 



jfc=i 



where Ug is some neighbourhood of the point , is such that 



oo 




for all (t, x) G U, 



00 • 



k=l 



Since solutions of the Cauchy problem depends holomorphically on initial data, we see 
that the function F is holomorphic on the neighbourhood Uqq. The Jacobi matrix at the 

8 



V.N.Gorbuzov 



Integral equivalence of multidimensional differential systems 



point (tgjXQ) with respect to x is the identity matrix, i.e., d^F{t,x)\ = E. Therefore 

Jacobian detd^F{t,x) for all {t,x) G Uqq . 

This implies that the coordinate functions F^: Uqo — M, z = 1, . . . , n, of the vector 
function F are functionally independent of x on the neighbourhood Uqq. 

By Lemma 1.1, 

F{t, x{t; {tQ , X J)) = x(io; {t, x{t; {Iq , xj))) = x^ . 

Hence the function F is a constant vector along a solution of system (TD). Then, by Theo- 
rem 1.1, the scalar functions 

F^: {t,x) ^ Fi{t,x) for all (t, x) G Uqq , i = 1, . . . ,n, (1.13) 

are first integrals on the neighbourhood Uqq of system (TD). ■ 

Lemma 1.3. Suppose the completely solvable system (TD) with Xg C°°{Il) has n fun- 
ctionally independent on a neighbourhood Uqq of the point {Iq^Xq) G 11 first integrals (1.13). 
Then for any first integral Uqq M of system (TD), we have ^'(i,x) = ^{F{t,x)) for all 
{t, x) G Uqq , where ^ is some scalar holomorphic function on the range EF of the vector 
function F: {t,x) — > [Fi{t, x), . . . , F^{t, x)) for all {t,x) G Uqq. 

Proof. The first integrals (1.13) of system (TD) has the form F-: {t,x) Xj^{tQ; {t, x)) 
for all {t, x) G Uqq , i = 1, . . . ,n. Jacobian det F{t, x) ^ for all {t, x) G Uqq . Therefore 
the function F at fixed t has the inverse function S and 

F{t,Sit,x)) = X for all (t, x) G Uqq . (1.14) 

The function $: {t,x) — > ^{t,S{t,x)) for all {t,x) G Uqq with the first integrals (1.13) 
is connected by the identity 

^(t, x) = $(t, F{t, x)) for all (t, x) G Uqq . 
Let us prove that the function $ is independent of t on the neighbourhood Uqq : 

d^,^{t,x) = for all (t,x)GUoQ, j = l,...,m. 
Differentiating the identity (1.14) with respect to t, we get 

dt,F{t,S{t,x)) + d^F{t,S{t,x))dt,S{t,x) = for all (i,x)GUQQ, j = l,...,m. 
The functions (1.13) are first integrals of system (TD). Hence, 

di.,F(t, x)= -d^ F(t, x) X^ (t, x) for all {t, x) G Uqq , j = l,...,m, 

where the vector functions X'' : (t,x)— > (Xy(t,x), . . . ,X^j{t,x)) for all (t,x)Gn, j = 1, . . . ,m. 
Then 

F{t, S{t, x)) {dt, S{t, x) - X^ {t, S{t, x))) = for all (t, x) G Uqq , j = l,...,m. 
Therefore, 

S{t, x) = X'' (t, S{t, x)) for all (t, x)gUqq, j = l,...,m, 

since F is a nonsingular matrix on Uqq . 

Taking into account the function ^' is a first integral of system (TD), we obtain 

dt.^t, x) = dt,^{t, Sit, x)) + d, *(i, Sit, x)) dt,Sit, x) = 

= at.^(i,5(t,x)) + 5^^(i,5(t,x))X^(t,5(i,x)) =0 for all (t,x) G Uqq , j = l,...,m. ■ 
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Prom lemmas 1.2 and 1.3, we have the following 

Theorem 1.3. The completely solvable system (TD) with X e C^{Il) on a neighbour- 
hood of any point of the domain H has a basis of first integrals of dimension n. 

Example 1.3. The completely solvable autonomous total differential system 

da;i=dti, = dt2, dx^ = d^^g{x^,X2)dt^ + d^_^g{x-^,X2)dt2, (1.15) 

where the scalar function g is holomorphic on a domain D C M^, has the basis of first 
integrals on the simply connected domain 11' = X D X M 

: (i, x) ^ t]^- x^ for all (t, x) G 11', F2 : (t, x) ^ ^2 ~ ^2 fo^' ^-H (t, x) G W, 

(1.16) 

F3 : {t, x) —>■ g{xi,X2) — x^ for all {t, x) G 11 . 

Indeed, by Definition 1.1, the functions (1.16) are first integrals on the domain 11' of 
system (1.15). The Jacobi matrix J(F^, F2, F3 ; t, .x) has rankJ{Fi,F2,F^;t,x) = 3 for all 
{t,x) G n'. Thus the first integrals (1.16) are functionally independent on the domain H'. 

Therefore, by Theorem 1.3, the first integrals (1.16) are an integral basis on the domain 
n' of the completely solvable autonomous total differential system (1.15). 



2. First integrals for linear homogeneous system 
of partial differential equations 

2.1. Basis of first integrals 

Definition 2.1. A scalar function F G C^{G') is said to be a first integral on a domain 
G' CG of system (d) with 2,^ G C(G), j = l,...,m, if 

StjF{x) = for all X EG', j = l,...,m. (2.1) 

Using k scalar functions 

Fg-. X ^ Fg{x) for all X G G', s = l,...,k, 

(2.2) 

GCi(G'), s = l,...,k, G'cG, 

we form the vector function 

F: x^ {F^{x),...,F,^{x)) for all x G G' (2.3) 

with range EF cR'^. 

Theorem 2.1. Let the functions (2.2) be first integrals on a domain G' G G of system 
{d) with 2,j G C{G), j = 1, . . . ,m. Then the function 

*: x^$(Fi(x),...,Ffc(x)) forallxeG', (2.4) 

where arbitrary function $ G C^(EF), is a first integral on the domain G' of system {d). 
Proof. By Definition 2.1, we have 

^^.^^(x) = for all X G G', j = 1, . . . ,m, s = 1, . . . ,k. 

Then 

k 

2.^F{x)) = Y,dp^F)\p^p^^^&jF,ix) = forallxGG', j = l,...,m. 

s=l 

Therefore the composite scalar function (2.4) is a first integral of system (d). ■ 
This theorem expresses functional ambiguity of first integrals for a linear homogeneous 
system of partial differential equations (see Subsection 1.2). 
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Definition 2.2. A set of functionally independent first integrals on a domain G' C G 
of system {d) with £,j G C{G), j = 1, . . . ,m, is called a basis of first integrals {integral 
basis) on the domain G' of system (d) if for any first integral ^ on the domain G' of 
system (d), we have ^{x) = ^{F{x)) for all x € G', where $ is some function of class 
C^(EF), EF is the range of function (2.3). The number k is said to be the dimension of 
basis of first integrals on the domain G' of system (d). 

Suppose m = n. Then the coefficient matrix 

u{x) = \\u..{x)\\ forallxGG (2.5) 

of system (d) is a square matrix of order n. Since the operators (0.1) are not linearly bound 
on the domain G, we see that the matrix (2.5) nearly everywhere on the domain G is 
nonsingular. Let G' C G be a domain such that rank«(x) = n for all x e G'. Then the 
system (d) on the domain G' is equivalent to the system of n differential equations 

dxiV = 0, i = l,...,n. 

The function y: x ^ C for all x G G', where C is arbitrary real constant, is a first 
integral on the domain G' of this system. 
Thus we obtain the following 

Property 2.1. If m = n, then all first integrals of system (d) are identical constants. 
Therefore the basic object of our research is the system (d) with m < n. 

2.2. Incomplete system 

Suppose the coordinate functions ti^j , j = 1, . . . ,m, i = 1, . . . ,n, of the linear differential 
operators (0.1) are a sufficient number of times continuously differentiable or holomorphic on 
the domain G, i.e., £j e C*^(G) or 2j G C~(G), j = l,...,m, where keN. 

Lemma 2.1. // the function y G C'^{G') is a first integral on a domain G' d G of the 
system of equations 

£i(x)y = 0, £2(^)y = 0, 

then this function is a first integral on the domain G' of the linear homogeneous first-order 
partial differential equation 

[£i(x),£2(2;)]y = 0. 

Proof [92, pp. 520 - 521]. The operations of Poisson brackets and commutators [91, p. 
165] on twice continuously differentiable functions coincide. This implies that 

[£i , £2] y{x) = £1^2 y{x) - £2-21 y{x) = O for all X G G'. ■ 

The following lemma is an immediate consequence of Lemma 2.1. 

Lemma 2.2. // the function F G C'^{G') is a first integral on a domain G' C G of 
system (d) with 2j G G^(G), j = 1, . . . ,m, then this function is a first integral on the domain 
G' of the linear homogeneous system of partial differential equations 

itj{x)y = Q, j = l,...,m, [&j^{x),Zi {x)]y = 0, 

(2.6) 

u = l,...,mi, H = l,...,m2, ^ m, m2 ^ m, , li^ e {1, . . . ,m}. 

Let the system (d) be complete. If to add to the system (d) at least one partial differ- 
ential equation of the form 

[£^.^(x),£;^(x)]y = 0, i,, /^G{l,...,m}, (2.7) 

then tlie system of partial differential equations (2.6) is built on the basis of the linearly bound 
on the domain G operators. 
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If the system (d) is incomplete, then some operators 

£^e{x)= T,ee{l,...,m}, 

are not a hnear combination of the operators £j , j = l,...,m. Adding the equation 
£^g(x) y= to the system {d), we obtain the Hnear homogeneous partial differential system 

£,g{x)y = 0, s = l,...,ki, m<k^<n, (2.8) 

where the operators s = 1, . . . ,ki, are not linearly bound on the domain G. 
It is important to underline, that 

1) ki > m, i.e., to system (d) one equation is added at least; 

2) the supplement of system (d) to the system (2.8) is made by adding of the equations 
of the form (2.7); 

3) the systems (d) and (2.8) are integrally equivalent on a domain G' C G, i.e., each first 
integral on the domain G' of system (d) is a first integral on the domain G' of system (2.8) 
and on the contrary each first integral on the domain G' of system (2.8) is a first integral on 
the domain G' of system (d). 

If the system (2.8) is complete, then the process is completed. 

If the system (2.8) is incomplete, then we do the similar procedure to the system (2.8) 
and obtain one more system. 

Note that on each step of this procedure tlie number of equations of an incomplete system 
increases at least by one. This argument shows that after a finite number of steps we get 
either the complete system or the system of n equations. 

Note also such conditions. 

First note that the system (d) with m = n is a complete system. Indeed, this follows 
from that the set of n not linearly bound on the domain G from the n-dimensional space R" 
first order linear differential operators of n variables is a basis of linear differential operators 
on the domain G. 

Therefore we may say that any incomplete system can be reduced to the complete system 
by a finite number of steps of the described procedure. 

Secondly note that the incomplete system {&) is reduced to the integrally equivalent 
complete system by adding the linear homogeneous partial differential equations of the forms 

[ii^.^(x),ii;^(x)]y = 0, [£,^(x), [£^.^(x),£,^(x)]] y = 0, 
£^^(x), \^^^{x), [St.^{x),Sli^{x)\]]y = 0, . . . , 



(2.9) 

i/=l, ...,mi, fi = l,...,m2, ^ = 1, ...,m3, C = 1, ■ ■ ■ , "^4 > ■ ■ ■ > 
rrig^m, s = 1, 2, . . . , {1, . . . ,m} 3 , l^, a^, (3^, .... 
Using these notations, we can state the following 

Theorem 2.2. Any incomplete system (d) with £,j € C^{G), j = 1, . . . ,m, A; G N, can 

be reduced to the integrally equivalent on some domain G' G G complete system by adding 
the equations of the form (2.9) to the system (d). 
At this point, we may give 

Definition 2.3. We'll say that a number S is the defect of the incomplete system (d) 
if this system can be reduced to the integrally equivalent on some domain G' C G complete 
system by adding S equations of the form (2.9). 

In this definition we mean that the corresponding complete system to the incomplete 
system {d) is constructed on the base of not linearly bound on the domain G operators. 

It is obvious that the incomplete system (d) has the defect S such that < 6 ^ n — m. 
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We may assume that a complete system has the defect 5 = 0. Then any system (d) has 
the defect 6 such that ^ S ^ n — m. 

Example 2.1. The normal linear homogeneous system of partial differential equations 

2,{x)y = 0, £2(^)^ = 0, (2.10) 

with the linear differential operators 

Sii{x) = d^^ + Xg d^^ — X4 d^^ for all x G M^, 

i22(x) = d^^ + 2x3X5 d^^ + 2x4X5 d^^ + (1 — X3 — X4 + X5) d^^ for all x G 

is incomplete as the Poisson bracket 

£21 (x) = [£2(2;), £i(^)] = 2x3X4 d^^ + (1 - xi + xl - xi) d^^ + 2x4X5 d^^ for all x G 

is not the null operator. 

The system (2.10) is integrally equivalent on some domain G' C to the linear homo- 
geneous system of partial differential equations 

iii(x)y = 0, S.^ix)y = 0, 2^^{x)y = 0. (2.11) 

The Poisson bracket 

^i-M^) = ['Ci(x), £21(2^)] = 2x3X5 d,^ + 2x4X5 d,^ + (1 - xi - x| + xi) d,^ for all x G 

is not a linear combination of the operators ii^, £,2, &2i- Therefore the partial differential 
system (2.11) is incomplete. 

The system (2.10) is integrally equivalent on some domain G' C to the linear homo- 
geneous system of partial differential equations 

£i(.T)y = 0, £2(^)^ = 0, £21(^)2/ = 0, £i.2i(x)y = 0. (2.12) 

The Poisson brackets 

[£i(x),£i.2i(x)] = -£2i(x) for all x G M^ 

[£2i(x),£i.2i(x)] = 4x5 - 4x4 = 

" 1 _ ^2 _ "1.2 _ ^2 -^21(3^) - 1 ™2 _ 12 _ ^2 ^l;2l{x) all X G G, 

J- 5 3 ^4 5 

[£2(x),£2i(x)] = [5,^ +£i;2i(x),£2i(x)] = - [£2i(x),£i;2i(x)] for all X G 
[£2(x),£i.2i(x)] = [5,^+£i;2i(x),£i;2i(x)] =0 for all x G M^ 

where G is any domain from the set D = {x: x| -|- X4 -)- x| 7^ 1 } of the space M^. 

Thus the Poisson brackets of the operators £1, £2, £21, £1-21 are the linear combinations 
of these operators on the domain G. Therefore the system (2.12) is complete on the domain G. 

The system (2.12) is obtained by adding of two equations to the system (2.10). The 
incomplete system (2.10) is integrally equivalent on some domain G' from the set D to the 
complete system (2.12). Thus the incomplete system (2.10) has the defect 6 = 2. 

2.3. Complete system 

Property 2.2. A complete linear homogeneous system of partial differential equations is 
invariant under a holomorphism. 
Proof. Let the map 
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x:^^ ifii) for all ^ G C 



(2.13) 



be a holomorphism between the domain Q. and the domain G of the space W^. 
The expression £j(x)y(x) is invariant under the holomorphism (2.13): 

ii/x) = ^(0 for all ^e^, for all x G G, j = 1, . . . , m, (2.14) 

where z{^) = y{^{^)) for all ^ € O,. Using the transformation (2.13), we have the system 
(d) is reduced to the system 



£,(0^ = 0, j = l,...,m, 



(2.15) 



with the not linearly bound on the domain Q linear differential operators of first order 
£,j , j = 1, . . . ,m (because a holomorphism is bijective). 

Let us prove that the system (2.15) is complete on the domain fl under the condition the 
system (d) is complete on the domain G. 

Taking into account the identities (2.14), we have 



StAx)Zi{x)y{x)\ 



for all ^ G for all X G G, j = 1, . . . ,m, 1 = 1,..., m, 
where f;(x) = y{x) for all x E G, 1 = 1,..., m. Hence, 



2Ax),£i{x) y{x) 



x=<p{^) 



^jiO^^liO ^(0 for all C e for all xEG, 



(2.16) 



j = l,...,m, l = l,...,m. 
Since the system (d) is complete, we see that the identities (0.12) are fulfilled. Therefore, 



E V(x)£,(x)y(x) 



x=(p{(,) 



Yl ^uiO ziO for alU G for all X G G, j = l,...,m, 1=1,..., 



m. 



Property 2.3. A jacobian linear homogeneous system of partial differential equations is 
invariant under a holomorphism. 

The proof is analogous to the proof of Property 2.2 if we take into account that the 
identities (0.14) for the jacobian system (d) are valid. ■ 

Property 2.4. The complete system (d) in a neighbourhood of any point x G G that 
satisfies dct[|V-'j/(x) [| ^ can be reduced to an integrally equivalent on some domain G' C G 
complete system by the nonsingular on the domain G linear transformation of the operators 



£j(x) = ^ i^jiix) %{x) for all X G G, j = l,... 



rn, 



(2.17) 



1=1 



where the linear differential operators of first order Otp I = 1, . . . ,m, and the scalar functions 
tpji : G — > M, j = 1, . . . ,m, I = 1, . . . ,m, are holomorphic on the domain G. 

Proof. The linear transformation (2.17) is nonsingular. Hence the linear differential oper- 
ators I = 1, . . . ,m, can be presented as the linear combinations of operators (0.1): 



iyiiix) = Y eij{x)£,j(x) forallxGGcG, l = l,...,m. 



(2.18) 
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and from the system (d), we obtain the system 

m 

^ VjjW5t/Wy = 0, j = l,...,m. (2.19) 

1=1 

The domain G in the expression (2.18) such that det||^^.;(x)|| for all x e G. 
The system (2.19) disintegrates on the system of equations 

%{x)y = 0, l = l,...,m, (2.20) 

where the linear differential operators of first order 9Tp I = 1, . . . ,m, are not linearly bound 
on the domain G. 

Prom the notion (2.19) and that the matrix ||'0^;(a;)|| of order m is nonsingular on the 

domain G C G it follows that the system (2.20) is integrally equivalent on some domain 
G' CG to the system (d). 

Let us show that the system (2.20) is complete. 

Using the operator identities 

[fix) £ .(x), g{x) £;(x)] = fix) g{x) [£ .(x), £;(x)] + /(x) £.gix) £,(x) - g{x) ilj{x) £ .(x) 
forallxeG, f e C^{G), geC^{G), j = l,...,m, 1 = 1,. ..,m, 

and 

[2j{x)+&i{x),&^{x)] = [2j{x),&^{x)] + [ii;(x),£^(x)] for all xEG, 
j = l,...,m, l = l,...,m, C = l,...,'m, 
and the representations (2.18), we get 

mm m 

[^^(x),^,(x)] = Aj,^,(x)[£,(x),£K^)] + j;i?.^.(x)£,(x) 

j=l 1=1 s=l 

for all X G G, n = 1, . . . ,m, v = 1, . . . ,m. 

Prom here using the decompositions (0.12) and the transformation (2.17), we obtain the 
Poisson brackets [9T^(.t), 9T^(.t)] for all .x G G, fi = l,...,m, v = l,...,m, are linear 

combinations on the domain G of the operators , j = 1, . . . ,m. This means that the 
system (2.20) is complete. ■ 

Prom Property 2.4, we have the following 

Property 2.5. // the system (d) is complete, then the system 

^j{x)y = 0, j = l,...,m, (2.21) 
where the linear differential operators of first order 

m 

2)^(x) = ^ Vji{x) 2,i{x) for all x G G, j = 1, . . . , m, (2.22) 

1=1 

the functional square matrix v{x) = \\'Vji{x)\\ of order m is nonsingular on the domain G, 
is also complete and integrally equivalent to the system (d) on a neighbourhood of any point 
X G G that satisfi.es dctf(x) ^ 0. 

Theorem 2.3. The complete system (d) can be reduced to an integrally equivalent on 
some domain G' G G complete normal system by the nonsingular on the domain G linear 
transformation of operators (0.1) {under this transformation we may have an restriction of 
the domain G). 
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Proof. Let the system (d) be complete. Then the square matrix u{x) = \\uj^{x)\\ for all 
X G G of order m (this matrix is obtain from the mxn matrix u{x) = \\uj^{x)\\ for all xGG 
by taking out the first m columns) is nonsingular on the domain G (since rank it(a;) = m 
almost everywhere on G, we see that it always can be received by renumbering variables). 
Therefore there exists a linear nonsingular transformation of operators (0.1) such that the 
system (d) can be reduced to the normal system (Nd). 

By Property 2.5, this normal system is complete. I 

Prom Theorem 2.2, we get the following 

Theorem 2.4. The complete system (d) can be reduced to an integrally equivalent on some 
domain G' <ZG jacobian system by the nonsingular on the domain G linear transformation of 
operators (0 .1) {under this transformation we may have an restriction of the domain G). 

Using Property 2.4 (or Property 2.5) and the process of building of complete normal system 
(see the proof of Theorem 2.3), we get the following theorem about integral equivalence of 
the complete system (d) and a complete normal system. 

Theorem 2.5. Suppose the complete system (d) has the nonsingular in the domain G 
square matrix u of order m (this matrix is obtain from the mxn matrix u{x) = \\uj^{x)\^\ 
for all X & G by taking out the first m columns). Then the complete system {d) can be 
reduced to the complete normal system (Nd) and these systems in a neighbourhood of any 
point X E G that satisfies det u{x) ^ are integrally equivalent. 

Theorems 2.2 and 2.5 are the substantiation of the following notion [1, p. 48]. 

Definition 2.4. We'll say that a domain H C G is a normalization domain of system 
(d) if the system (d) in a neighbourhood of any point of the domain H can be reduced to an 
integrally equivalent complete normal system. 

Notice that a normalization domain of an incomplete system is a normalization domain 
of an integrally equivalent complete system (see Theorem 2.2). 

In general case a normalization domain is ambiguous. This normalization domain depends 
on zeroes of determinants of square matrices of order m (these matrices are obtain from the 
matrix u{x) = ||^ji(^)||^xn ^^'^ x E G by taking out m columns). 

Example 2.2. Consider the linear homogeneous system of partial differential equations 

£i(x)y = 0, £^{x)y = 0, (2.23) 
where the linear differential operators of the first order 

S,i{x) = X]^ d^^ + X2 d^^ + x^ d^^ + x^ d^^ + x^ d^^ for all x G R^, 

£2(2;) = Xi d^^ + X2 d^^ + X3 d^^ + x\ d^^ + X5 d^^ for all x G M^. 
The Poisson bracket 

[£i(a;), £2(x)] = xl d^^ + d^^ = Q^^i^) for all xER^ 

is not a linear combination of the operators £,1 and £2- Hence the system (2.23) is incomplete. 
Using the operator £^2) g^t the system (2.23) is reduced to the system 

£i(,x)y = 0, £2(x)y = 0, £,2ix)y = 0. (2.24) 

Since the Poisson brackets 

[£i(x), £12(3;)] = £12(2:) forallxGM^ [£2(x), £i2(a;)] = D for all x G R^ 

we see that the system (2.24) is complete. 

Therefore the incomplete system (2.23) has the defect 5 = 1. 

Prom the second equation of system (2.24) by virtue of the third equation of this system 
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it follows that 

Xid^^y + X2d^^y + x^d^^^y = 0. (2.25) 
Then from the first equation of system (2.23), we have 

^A^y + H^x^y = 0- 

Prom this equation and the third equation of system (2.24), we obtain the equalities 

dx^y = 0, d^^y = 0. 

By solving the equation (2.25) for we reduce the system (2.24) to the normal system 

dx^y= - ^2^1^ dx^y - d^^y^ 9^^y = o, d^^y = o- (2-26) 

The complete normal system (2.26) is integrally equivalent to the complete system (2.24) 
and to the incomplete system (2.23) on a normalization domain C {x: 0} C M^. 

It is readily seen that the systems (2.23) and (2.24) have else two normal forms, which 
can be obtain by solving the equation (2.25) for d^^y and d^^y- 

dx^y= -XiX2^ d^^y -x^X2^ d^^y, d^^y = 0, d^^y = (2.27) 

and 

dx,y= -xix^^d^^y-x^x^^d^^y, d^^y = 0, d^^y = 0. (2.28) 

The complete normal system (2.27) is integrally equivalent to the complete system (2.24) 
and to the incomplete system (2.23) on a normalization domain H2 C {x: X2 7^ 0} C M^. 

The complete normal system (2.28) is integrally equivalent to the complete system (2.24) 
and to the incomplete system (2.23) on a normalization domain C {x: x^ 0} C M^. 

2.4. Dimension of integral basis 

The system of total differential equations 

m 

dXg = — "Ujglx) dxj , s = m + 1, . . . ,n, (2.29) 

is associated to the normal linear homogeneous system of partial differential equations (Nd). 
For the system (2.29), as well as for the system (TD), the linear differential operators Xj , 
j = 1, . . . ,m, have the form 

Xjix) = &j{x) = d^, - m^{x) for all X G G, j = l,...,m, (2.30) 

where the operators 971^- , j = 1, . . . ,m, are given by (0.12). 

Definition 2.5. We'll say that a total differential system and a linear homogeneous system 
of partial differential equations are integrally equivalent on some domain if on this domain 
each first integral of the first system is a first integral of the second system and on the contrary 
each fi,rst integral of the second system is a first integral of the first system. 

If a total differential system and a linear homogeneous system of partial differential equa- 
tions are integrally equivalent on the domain 11', then these systems have the same integral 
basis on this domain (see Definitions 1.2, 2.2, and 2.5). 

Prom the connections (2.30) it follows that the identities (1.2) for the function F: G' — > M 
coincide with the identities (2.1) for the function F: G' — > M. By Definitions 1.1 and 2.1, we 
get the next theorem about the integral equivalence of a linear homogeneous system of partial 
differential equations with a system of total differential equations. 
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Theorem 2.6. The scalar function F: G' is a first integral on the domain G' d G 

of the normal linear homogeneous system of partial differential equations if and only if this 
function is a first integral on the domain G of the total differential system (2.29). 

Using the definitions of the complete and jacobian hnear homogeneous systems of partial 
differential equations {d) and the definition of the completely solvable total differential system 
(TD), we get the relation between these notions. 

Theorem 2.7. The normal linear homogeneous system of partial differential equations (Nd) 
is complete (jacobian) if and only if the total differential system (2.29) is completely solvable. 

By Theorem 2.6, using the definitions of an integral basis for a total differential system 
(Definition 1.2) and for a linear homogeneous systems of partial differential equations (Defi- 
nition 2.2), we have the relation between integral bases of the systems (Nd) and (2.29). 

Theorem 2.8. The set of functions (2.2) is a basis of first integrals on a domain G' (Z G 
of the normal linear homogeneous system of partial differential equations (Nd) if and only if 
this set of functions (2.2) is a basis of first integrals on the domain G' of the system of total 
differential equations (2.29). 

Prom Theorem 1.3 it follows that the dimension of an integral basis for the completely 
solvable total differential system (2.29) is equal n — m. Then, using Theorems 2.7 and 2.8, we 
obtain the dimension of an integral basis for the complete normal linear homogeneous system 
of partial differential equations (Nd). 

Theorem 2.9. The complete (jacobian) normal linear homogeneous system of partial 
differential equations (Nd) with Tlj G C°°{G) on a neighbourhood of any point from the 
domain G has an integral basis of dimension n — m. 

Example 2.3. The complete (jacobian) normal linear homogeneous system of partial 
differential equations 

\y= -9a;^y-d^,9{xi,X2)d^^y, d^^y = - d^^y - d^^g(x^ ^x^) d^^y, (2.31) 

where the scalar function g G C°°(D), D C M^, is associated to the completely solvable 
system of total differential equations (1.15). 

The system (2.31) has the form (Nd) with rn = 2, n = 5. 

In Example 1.3 an integral basis on the domain H' = x D x M for system (1.15) was 
constructed. This integral basis is three functionally independent on the domain 11' first 
integrals (1.16) of system (1.15). 

By Theorem 2.6, the functions (1.16) are first integrals on the domain 11' of the partial 
differential system (2.31). 

Prom Theorems 2.8 and 2.9 it follows that the complete normal system (2.31) has an 
integral basis on the domain 11' of dimension n — m = 5 — 2 = 3. This integral basis is the 
functionally independent on the domain 11' first integrals (1.16). 

Using the definition of the normalization domain (Definition 2.4), Theorems 2.5 (to the 
effect that the complete system (d) can be reduced to an integrally equivalent complete 
normal system) and 2.9 (about the dimension of an integral basis of the complete normal 
system (Nd)), we obtain the dimension of an integral basis for the complete system (d). 

Theorem 2.10. The complete linear homogeneous system of partial differential equations 
(d) with £j G C°°(G), j = I, . . . ,m, on a neighbourhood of any point of a normalization 
domain has a basis of first integrals of dimension n — m. 

Example 2.4. We consider the linear homogeneous system of partial differential equations 

£i(x)y = 0, £2(^)^ = (2-32) 
with the linear differential operators of first order 

iJ^(x) = d^^ + d^^ + d^^ for all x G M^, £2(2^) = + ^2 + "^2, ^x^ for all x G M^. 
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Since the Poisson bracket 

[&^{x),£2{x)] = S,^{x) for all x G M^, 

we see that the system (2.32) is complete. 

By Theorem 2.10, the system (2.32) has an integral basis of dimension n — m = 3 — 2 = 1. 

Using the definition of a first integral (Definition 2.1), we obtain a first integral of system 
(2.32) and consequently an integral basis of system (2.32). 

Thus the function 

Fi2- X ^ {x2- x^){x^- x^)''^ for all x G ' = {x: X2 x^}, 

is a basis of first integrals for the system (2.32) on any domain from the set Di2- 
Similarly the function 

F^^ : X (.T3 — 3^2) (-^1 ~ ^3)^^ fo'^ ^ ^ -^13 ' ^13 — ^3 7^ ^1}' 

is a basis of first integrals for the system (2.32) on any domain from the set Di^. 
In the same way the function 

^23:0;^ {x^ - a;i)(a;2 - x^)'^ for all xe D^^, -D23 = {x- x^^ X2}, 

is a basis of first integrals for the system (2.32) on any domain G23 from the set 1)23 • 

Each of the integral bases , -^13, and F23 of the linear homogeneous partial differential 
system (2.32) defines the family of integral surfaces of this system which consists of the planes 
C^Xi + C2X2 — {Ci + C2)x3 = 0, where and C2 are arbitrary real constants. 

Using the notion of defect for a system (Definition 2.3), the theorem of the reduction of 
an incomplete system to an integrally equivalent complete system (Theorem 2.2), and the 
theorem of dimension for an integral basis of a complete system (Theorem 2.10), we get the 
theorem of dimension for an integral basis of an incomplete system. 

Theorem 2.11. Suppose the incomplete linear homogeneous system of partial differential 
equations (d) with 2,j G C°°(G), j = 1, . . . ,m, has the defect 6. Then this system on a neigh- 
bourhood of any point of a normalization domain has an integral basis of dimension n—m—S. 

Using Property 2.1 and the procedure of the reduction of an incomplete system to a 
complete system (see Subsection 2.2), we obtain 

Corollary 2.1. Suppose the incomplete linear homogeneous system of partial differential 
equations {d) has n — 1 equations with n unknowns. Then first integrals of this system are 
only arbitrary constants. 

Recall that the complete system (5) has the defect 5 = 0. Using this notation. Theorems 
2.10 and 2.11, we can state the following 

Theorem 2.12. Suppose the linear homogeneous system of partial differential equations 
(d) with £,j G C°°{G), j = 1, . . . ,m, has the defect 6, ^ S ^ n — m. Then this system 
on a neighbourhood of any point of a normalization domain has a basis of first integrals of 
dimension n — m — S. 

Prom Theorem 2.12, we obtain the completeness criterion for linear homogeneous system 
of partial differential equations. 

Theorem 2.13. The linear homogeneous system of partial differential equations (d) with 
&j G C°°{G), j = 1, . . . ,m, is complete if and only if this system on a neighbourhood of any 
point of a normalization domain has a basis of first integrals of dimension n — m. 

The following agreement is needed for the sequel. 

Agreement 2.1. By ^^^{x) for all x E G, v = l,...,p, denote linear differential 

operators of first order, which constructed on the base of the operators (0.1), such that the 
operators , . . . , £^ , £^ , . . . ,£,.^ are not linearly bound on the domain G and the equations 
£,l{x)y = 0, v = 1, . . . ,p, have the forms (2.9). 
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Along with the system (d) we'll consider the linear homogeneous system of partial dif- 
ferential equations 

£j{x)y = 0, j = l,...,m, £*(x)y = 0, u = l,...,p. (2.33) 

The system (2.33) is constructed on the base of system (d) according to Agreement 2.1. 

By Lemma 2.2 and in accordance with Agreement 2.1, we obviously have 

Theorem 2.14. The system {d) and the system (2.33), which constructed on the base of 
system (d) according to Agreement 2.1, are integrally equivalent on some domain G' C G. 

Prom Theorem 2.14, we get the following assertion for an integral basis. 

Corollary 2.2. A set of scalar functions is an integral basis of dimension r on a domain 
G' d G for the partial differential system (d) if and only if this set is an integral basis of 
dimension r on the domain G' for the partial differential system (2.33), which constructed 
on the base of system (d) according to Agreement 2.1. 

The system (2.33) under the condition p = is the system (d). 

If p = 6, where 5 is a defect of system {d), then the system (2.33) is complete. 

Using the pocedure of the reduction of an incomplete system to a complete system (see 
Subsection 2.2), Agreement 2.1, Theorem 2.10, and Corollary 2.2, we clearly have 

Theorem 2.15. The scalar functions : G' ^ R, t = I, . . . ,m — n — 5, are a basis of 
first integrals on the domain G' <Z G for system (d) with £j G C°°{G), j = 1, . . . ,m, and 
with the defect 6, ^ 6 ^ n — m, if and only if this functions are a basis of first integrals on 
the domain G' for the complete system (2.33), where p = S. 

Example 2.5. In accordance with the definition of a first integral (Definition 2.1), we 
obtain the scalar function 

F:x^ x^il + xl + xl + xl)-'^ for all x G (2.34) 

is a first integral of system (2.10). 

In Example 2.1 we proved that the system (2.10) is incomplete and has the defect 5 = 2. 
By Theorem 2.11, an integral basis of system (2.10) has the dimension n—m—5 = 5—2—2 = 1. 

Therefore the first integral (2.34) is an integral basis on the space of the incomplete 
system (2.10). 

By Corollary 2.2 and Theorem 2.15, it follows that the scalar function (2.34) is an integral 
basis on the space both the incomplete system (2.11) and the complete system (2.12). 

Example 2.6. In Example 2.2 it has been shown that the system (2.23) is incomplete 
and has the defect 6 = 1. 

By Theorem 2.12, the incomplete system (2.23) on a neighbourhood of any point of a 
normalization domain has a basis of first integrals of dimension n — m — 5 = 5 — 2 — 1 = 2. 

Prom the definition of a first integral (Definition 2.1) it follows that the complete normal 
system (2.26) has the first integrals 

: X X2Xi^ for all x & and F^^ : x x^x^^ for all x & H^, (2.35) 

where is any domain from the set {x: x^ ^ 0} of the space M^. 

By Theorem 2.15, the scalar functions (2.35) are a basis of first integrals on any domain 
C {x: x^ ^ 0} of the complete normal system (2.26), of the complete system (2.24), and 

of the incomplete systems (2.23). 
Similarly, the scalar functions 

F-^2' X — > XiX2^ for all x G H2 and F32 : x — ^ x^X2^ for all x G H2 (2.36) 

are a basis of first integrals on any domain H2 C {x : X2 ^ 0} of the complete normal system 
(2.27), of the complete system (2.24), and of the incomplete systems (2.23). 
The scalar functions 
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: X — > x^Xg for all x € and -F23 '■ x ^ ^^2^3 fo^' '^^^ ^ ^ -^3 (2.37) 

are a basis of first integrals on any domain C {x : Xg 7^ 0} of the complete normal system 
(2.28), of the complete system (2.24), and of the incomplete systems (2.23). 

Each of the integral bases (2.35), (2.36), (2.37) for the incomplete system (2.23) (the com- 
plete system (2.24)) defines the two families of the integral surfaces, which consistes of planes, 
for this system respectively 



P21 


= {x: 




+ C2X2 


= 0} 


and 


P31 


= {x: 




+ C4X1 


= 0}, 


P12 


= {x: 




+ C2X2 


= 0} 


and 


P32 


= {x: 


C3X3 


+ C4X2 


= 0}, 


Pl3 


= {x : 


CiXi 




= 0} 


and 


P23 


= {x: 


C2X2 


+ C4X3 


= 0}, 



where C^, . . . , C4 are arbitrary real constants. 

3. Dimension of integral basis for not completely solvable 
total differential system 

The normal on the domain 11 C R"*''"" linear homogeneous partial differential system 

n 

^t.y = -^^ij{i^x)d^,y, j = l,...,m, (3.1) 

i=l 

is associated to the system of total differential equations (TD). 
By Definitions 1.1, 2.1, and 2.5, we obtain 

Theorem 3.1. The total differential system (TD) is integrally equivalent on some domain 
n' C n to the normal linear homogeneous system of partial differential equations (3.1). 
Prom Theorems 2.7 and 3.1, we have 

Theorem 3.2. The total differential system (TD) is completely solvable if and only if the 
normal linear homogeneous system of partial differential equations {3 .1) is complete (jacobian). 
Using Theorems 1.3, 2.9, 3.1, and 3.2, we can state the following 

Theorem 3.3. // the total differential system (TD) with X G C°°(n) is completely 
solvable {the normal linear homogeneous system of partial differential equations (3.1) with 
Xj^j G C°°(n). i = 1, . . . , n, j = 1, . . . ,m, is jacobian), then the systems (TD) and (3.1) have 
the same integral basis of dimension n on some domain IT' C H. 

The total differential system (1.15) is completely solvable. The jacobian normal linear 
homogeneous system of partial differential equations (2.31) is associated to the system (1.15). 
The integral basis of these systems was built in Examples 1.3 and 2.3. 

Example 3.1. Let us consider the total differential system 

dxi = Xi{xi + l){dti + dt^)-, dx2 = X2{xi + 2){dti + dt2), 

(3.2) 

dxg = x^{xi + 3) dti + x^{xi + 5) dt2- 
The normal linear homogeneous system of partial differential equations 

Xi(t,x)y = 0, X2(t,x)y = 0, (3.3) 
where the linear differential operators of first order 

Xi{t, x) = + xi(xi + 1) d^^ + X2(xi + 2) d^^ + X3(xi + 3) d^^ for all {t,x) G 

X2{t, x) = + xi(xi + 1) d^^ + X2(xi + 2) d^^ + X3(xi + 5) d^^ for all {t, x) G R^ 

is associated to the total differential system (3.2). 
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Since the Poisson bracket [Xi{t,x),3c2(t,x)] = O for all {t,x) G M^, we see that the 
system (3.2) is completely solvable and the system (3.2) is complete. 

The systems (3.2) and (3.3) have the same integral basis of dimension r = 3. 

The functionally independent on any domain Hi C 'E^ = {{t, x): x-y^ 0} scalar functions 

: [t, x) x^ixi + 1) Xi^ for all {t, x) G , 
F2 : (t, x) (xi + 1) Xi^ exp(ii + ig) for all {t, x) € , (3.4) 
F3 : {t, x) ^ Xgx^^ exp( - 2{ti + 2*2)) for all {t, x) G 

are a basis of first integrals on the domain both for the system (3.2) and the system (3.3). 

By the assumption of functional ambiguity of first integrals (Theorems 1.2 and 2.1), we 
can build integral bases of the systems (3.2) and (3.3) which are different from the integral 
basis (3.4). Note also that domains of definition for these integral bases can contain the set 
of points {ti ,0,X2 ,x^). 

For example, the scalar functions F^^, F^^, and ^3"^ on any domain H2 from the set 
H2 = {(t, x) : Xi 7^—1, .T2 / 0, ^ 0} are a basis of first integrals on the domain H2 both 
for the system (3.2) and the system (3.3). 

If the system (TD) is not completely solvable, then the system (3.1) is incomplete. In this 
case, we reduce the system (3.1) to the complete system. We obtain the defect S, < 6 ^ n, 
and a normalization domain of system (3.1). By Theorem 2.12, we have the following 

Theorem 3.4. The not completely solvaMe total differential system (TD) with X€ C°°(n) 
on a neighbourhood of any point of a normalization domain for the linear homogeneous (in- 
complete) system of partial differential equations (3.1) has an integral basis of dimension n—S, 
where 6 is the defect of system (3.1). This basis of first integrals for system (TD) is also a 
basis of first integrals for system (3.1). 

Example 3.2. The normal linear homogeneous system of partial differential equations 

Xi{t,x)y = 0, X2{t,x)y = 0, (3.5) 

where and are the linear differential operators of first order (1.4) and (1.5) respectively, 
is associated to the total differential system (1.3). 
The Poisson bracket 

Xi2{tiX) = \Xi{t,x),X2(t,x)\ = tiX2X^d^^ + x^(l + Xiti — x\ti ) dj.^ + 

+ X2( - 2 - 2xiti^ + 2x1 + + ^3) for all (t x) G D 

is not the null operator. Therefore the system (1.3) is not completely solvable and the system 
(3.5) is incomplete. 

Since the Poisson bracket 

^l;12(^' ~ \^l{tiX)-,Xi2{tjX)\ = 
= - 2xSl - - Xj + Xi) d^^ + X2X2{1 - 5Xiti^ + x\ tX^) d^^ + 

+ (2 + 4xit^^ - 2x\ - 4xf tl"^ + 2x\ tl'^ + Xix\ t^^ + ?,x\x\ t^'^ + 2x\ - bxl - 

- 2x^x1 ti^ + 2x1x1 ti^ + 2x\xl - 2x1) ^x^ for all (*, x) G 

is not a linear combination of the operators X-y^ X2, X12, we see that the linear homogeneous 
system of partial differential equations 

Xi{t,x)y = 0, X2{t,x)y = 0, Xi2{t-ix)y = (3.6) 

is incomplete. 
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In Example 1.1, we proved that the function (1.6) is a first integral on a domain 11 C D 
of system (1.3). Therefore the incomplete system (3.6) has the defect 6 = 1 and the function 
(1.6) is an integral basis on a domain 11 C D of system (3.6). 

Thus the function (1.6) is an integral basis on a domain H c D both for the not completely 
solvable total differential system (1.3) and the incomplete normal system of partial differential 
equations (3.5) with the defect 6 = 2. 

Recall that a complete linear homogeneous system of partial differential equations has the 
defect 6 = 0. Hence for the completely solvable total differential system (TD) and for the 
not completely solvable total differential system (TD) we have the following assertion 

Theorem 3.5. The system (TD) with XGC°°(ir) and the associated system (3.1) to the 
system (TD) have the same integral basis of dimension n — 6 on a neighbourhood of any point 
of a normalization domain for system (3.1), where 6, ^ 6 ^ n, is the defect of system (3.1). 

Using this notation, we can state the definitions of defect and of normalization domain 
for completely solvable system (TD) and for not completely solvable system (TD). 

Definition 3.1. The total differential system (TD) has the defect 6, ^ 6 ^ n, where 6 
is the defect of the associated linear homogeneous system of partial differential equations (3.1). 
Thus a normalization domain of system, (3.1) is called a normalization domain of system (TD). 

If the system (TD) is completely solvable, then a normalization domain of this system is 
the domain of complete solvability for system (TD). 

By Definition 3.1 and Theorem 3.5, we obtain 

Theorem 3.6. Suppose the total differential system (TD) with X G C°°(n) has the defect 
6, ^ 6 ^ n. Then this system has an integral basis of dimension n — 6 on a neighbourhood 
of any point of a normalization domain. 

Example 3.3. The total differential system 

dxi = Xi dti + xf dt2 , dx2 = x\ dt^ + dt2 (3.7) 
induces the linear differential operators of first order 

{t, x) = + x^ d^^ + x\ d^^ for all (t, x) G M^, 

X2(t, x) = + xl d^^ + xl d^^ for all (t, x) G Ml 
Since the Poisson bracket 

dii2{t,x) = [Xi(t, x), X2(t, x)] =x\d^^ +^2'^x,^ for all (t, x) G 

is not the null operator, we see that the system (3.7) is not completely solvable. 
The linear differential operators of first order , X2 , ■> ^2- 12 > where 

^2; 12(^5 a;) = [X2(t,x),Xi2(t,a;)] = xld^^ for all (i,x) G M"^, 

are not linearly bound on R^. Therefore the associated incomplete normal linear homogeneous 
system of partial differential equations 

y.i{t,x)y = 0, X2{t,x)y = 

to the system (3.7) has the defect 6 = 2. 

Thus the system (3.7) has the defect 6 = 2. From n — 6 = 2 — 2 = it follows that the 
system (3.7) has no first integrals. 

Example 3.4. The normal linear homogeneous system of partial differential equations 

Xi{t,x)y = dt^y + xid^^y + {l + xi + 2x2)d^^y = 0, 

(3.8) 

X2{t, x)y = d^^y + 3xi d^^y + (x^ + 3x2) d^^y = 
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is associated to the not completely solvable total differential system (1.7). Therefore the 
system (3.8) is incomplete. 
The Poisson brackets 

Xi2(t, x) = [X^{t, x), X^it, x)] = (3 - Xi) for all {t, x) G M^ 
[Xi{t,x),Xi2{t,x)] = {xi — 6){3 — Xi)~^ Xi2{t,x) for all {t,x) E {{t,x): x-^ 3}, 

[X.2{t, x), Xi2{t, x)] = 9{xi — 3)~^ X.i2{t, x) for all {t, x) G {{t, x): x-^ 3}. 

Thus the system (3.8) has the defect 5 = 1. 

By Definition 3.1, the system (1.7) has the defect 6=1. By Theorem 3.6, the dimension 
of an integral basis of system (1.7) is n — 1 = 2 — 1 = 1. Therefore the function (1.8) is an 
integral basis on the space M*^ of system (1.7). 

By Theorem 3.5, the function (1.8) is an integral basis on the space of the system of 
partial differential equations (3.8). 

Example 3.5. The normal linear homogeneous system of partial differential equations 
(2.10) has the defect 5 = 2 and the integral basis (2.34). This system is associated to the 
total differential system 

dx^ — ^x^x ^ dx^ f dx ^ — dx-^ ~\~ ^x ^x ^ dx^ ^ 

(3.9) 

dx^ = — dxi + (1 — X3 — X4 + a:|) dx2 ■ 

Therefore the autonomous system (3.9) is not completely solvable and has the defect 
6 = 2. The function (2.34) is an autonomous first integral of system (3.9) and this function 
is an integral basis on the space of system (3.9). 

Suppose the system (TD) has the defect 6, ^ 6 < n. Then the associated system (3.1) 
to the system (TD) we reduced to the integrally equivalent complete system 

n 
i=l 

(3.10) 

n 

^ Xi^{t,x)d^,y = 0, u = l,...,6, 
1=1 

where the functions X*^ : 11 ^ R, i = 1, . . . ,n, u = 1, . . . ,6, are constructed on the base of 
the functions X^^, i = 1, . . . ,n, j = 1, . . . ,m, by the rule (2.9). 

We reduced the system (3.10) to a normal system and then for this normal system we 
build the associated total differential system 

m n 

dxk^ =^Gk^^{t,x)dtj + Gi,^k^{t,x)dxk^, 

j=l /i=n-(5+l 

(3.11) 

7 = l,...,n-5, k^, k^^ £{!,..., n}, k^^^k^, i = l,...,n, ^=l,...,n, i ^ ^. 

The system (3.11) is completely solvable on a normalization domain of system (3.10) (this 
normalization domain is a normalization domain both for the partial differential system (3.1) 
and the total differential system (TD)). 

The system (3.11) in relation to the system (TD) has the extended coordinate space Ot 
on 6 coordinates at the expense of <5 coordinates of the coordinate space Ox. 

Using the rearrangement of the dependent and independent variables in the system (TD), 
we get the system (3.11) has the form 
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m+S 

= ^ Hij{ti tm+s ,xi,..., Xn-s) dtj, i = l,...,n-d, (3.12) 

where = x^_s+^ , 1^ = 1,..., 5. 

Note also that the system (TD) and the system (3.12) are integrally equivalent (they 
have the same integral basis) on a normalization domain (the system (TD) on this domain is 
reduced to the system (3.12)) of system (TD). 

Theorem 3.7. The total differential system (TD) with the coordinate spaces Ot and 
Ox and with the defect S, ^ S < n, on a normalization domain is integrally equivalent to 
the completely solvable total differential system with the coordinate spaces Ot^ , • • • , t^^^g and 
Oxi , . . . , , where = •'^n-5+u ' = Ij • • • j ^ (accurate to numbering of the dependent 

variables x^ , . . . ,x^ in the system (TD)). 

Example 3.6. Adding the equation Xi2{t,x)y = (3 — Xi)d^^y = to the incomplete 
system (3.8), we get the integrally equivalent complete normal system 

The total differential equation 

dX]^ = Xi dti + 3xi dt2 + dx2 (3-14) 

is associated to the system (3.13). 

Therefore the not completely solvable total differential system (1.7) is integrally equivalent 
to the completely solvable total differential equation (3.14). Moreover, the system (1.7) and 
the equation (3.14) have the same integral basis, which is the first integral (1.8). 

In Example 1.2, we proved that the system (1.8) has no solutions. At the same time the 
integrally equivalent equation (3.14) to the system (1.8) has the general solution 

Xi : {ti , t2 , ^2) — > C ex.p{ti + 312) for all {t^ , t2 , X2) £ M^. 

Thus integrally equivalent total differential systems have the same first integrals. But this 
statement is not true for solutions of integrally equivalent systems. 

Example 3.7. The complete system of partial differential equations (2.12) is reduced to 
the normal systems: 

d v = d y = d _l-xl + xl + xl y = ^d v 

^x^y ^x^y ^x;^y 2x3X5 5 * x^ ^5^' 

d v = d y = d v= _'^-^l + ^l + ^5 q q y = ^Q y. 

^x^y ^> '-'x^y '-'x^y 2x3X4 ^4^' x^y^ 

dj. y = 0, dj. y = 0, y = - ^ d^y, d„ y = - ^ d„ y. 

^1 ^2 1-xl + xl + xl ^3 1-xl + xl + xl 

The completely solvable total differential equations 



1 — x| + X4 + x^ 



dxr, = dxi + dx2 1 dx^ + — dx^, (3.15) 

T3X5 Xg 



dX4 = dXi + dX2 - ^3 + 3^4 + 3^5 ^ + ^ (3 jg) 

2xo 



dxg = dxi + dx2 - _.^A,o dx^ - , _ 2 2 , ^2 (3.17) 

3 4 5 3 4 5 



are associated to these normal systems respectively. 
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Thus the not completely solvable total differential system (3.9) is: 

a) integrally equivalent on any domain C {x: ^ 0, ^ 0} to the completely 
solvable total differential equation (3.15) and they have the same integral basis on the domain 

, which is the first integral (2.34); 

b) integrally equivalent on any domain C {x: 2:3 7^ 0, 0:4 7^ 0} to the completely 
solvable total differential equation (3.16) and they have the same integral basis on the domain 
i?4 , which is the first integral (2.34); 

c) integrally equivalent on any domain G {x: 1 — x| + 2:4 + x| 7^ 0} to the completely 
solvable total differential equation (3.17) and they have the same integral basis on the domain 

, which is the first integral (2.34). 

4. First integrals for Pfaff system of equations 

4.1. Integrally equivalent PfafT systems of equations 

Definition 4.1. A scalar function F G C^{G') is said to be a first integral on a domain 
G' C G of system (Pf) with iVj G C{G), j = l....,m, if there exist the scalar functions 
ttj € C{G'), j = 1, . . . ,m, such that the total differential 

m 

dF{x) = aj{x)u}j{x) for all x G G'. (4.1) 

Let us introduce the equivalence relation on a set of Pfaff systems of equations. 
Definition 4.2. We'll say that two Pfaff systems of equations are integrally equivalent 

on some domain if on this domain each first integral of the first system is a first integral of 
the second system and on the contrary each first integral of the second system is a first integral 
of the first system. 

We claim that the demand that the 1-forms (0.2) are not linearly bound on the domain 
G is not narrow the set of all possible Pfaff systems of equations (Pf ) (from the point of view 
of the integral equivalence). Indeed, let the 1-forms (0.2) be linearly bound on the domain 
G. Then the functional matrix 

Hx) = |hji(a^)|Lxn fo'^ all X G G (4.2) 

has the rank 

rankiy(a;) = s{x), 1 ^ s{x) < min {m,n} for all x G G. 

Take s = min{s(x) : x G G} not linearly bound on a domain O C G 1-forms 

o;^. , G {1, . . . , m}, 1 = 1,..., s. (4.3) 

The domain ft is such that the complement ft on G has the null measure: /iCqH, = 0. 
Using the 1-forms (4.3), we get the new PfafF system of equations 

ijj^{x) = 0, G {l,...,m}, l = l,...,s. (4.4) 

Since the 1-forms (4.3) are not linearly bound and s = min{s(x): xGG}, we see that the 
PfafF systems of equations (Pf) and (4.4) have the same first integrals on the domain ft (by 
Definition 4.1). Thus the systems (Pf) and (4.4) are integrally equivalent on the domain O. 

Since CqQ has the null measure, we see that a class of systems (Pf) doesn't restrict. 

4.2. Integral basis 

Suppose the 1-forms (0.2) are not linearly bound on the domain G. Then the functional 
matrix (4.2) has ranku;(x) = m almost everywhere on the domain G. To be definite, assume 
that m ^ n. 
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If m = n, then, since the hnear differential forms (0.2) are not linearly bound on the 
domain G, it follows that the square matrix (4.2) of order n is nonsingular almost everywhere 
on the domain G. 

In this case the system (Pf ) on a domain fl C G by a nonsingular algebraic transformation 
can be reduced to the differential system 

dx^ = 0, i = 1, . . . ,n, 

where the domain J7 is such that /iC^Jl = 0. 

Whence, we obtain = Cj , i = 1, . . . , n, where C^^ , . . . , C„ are arbitrary real constants. 
Hence the functions 

F-: X ^ Xj^ for all x G J7, i = l,...,n, (4.5) 

are first integrals on the domain O of system (Pf). 
Thus the case m = n is singular and here we have 

Property 4.1. The system (Pf) with m = n has n functionally independent first 
integrals (4.5) on such a subdomain 0, of the domain G that hCqQ = 0. 

Theorem 4.1. Let the functions (2.2) be first integrals on a domain G' C G of system 
(Pf) with G G{G). Then the function (2.4) is also a first integral on the domain G' of 
system (Pf). 

Proof. By Definition 4.1, the functions (2.2) are first integrals on the domain G' of system 
(Pf) if and only if there exist the scalar functions a^j € C{G'), £, = 1, . . . ,k, j = 1, . . . ,m, 
such that the total differentials of functions (2.2) have the forms 

m 

dF^{x) = a^j{x)uj{x) for all x G G', = 1, • • • , ^- (4-6) 

Suppose $ is arbitrary scalar function from the space C^(EF), where F is the vector 
function (2.3). Then, using the identities (4.6), we get the total differential of function (2.4) is 



k 

d^{x) = d^{F,{x), . . .,F,ix)) = ^ dp^^iF, , . . • dF^ix) 

(=1 



= T.T.'^f^HF^,---,F,)^ a^^{x)u;^{x) for all x G G'. 

By this identity and Definition 4.1, it follows that the function (2.4) is a first integral on 
the domain G' of the PfafF system of equations (Pf). I 

It was shown in Theorem 4.1 that first integrals for a Pfaff system of equations are func- 
tional ambiguous. Thus the priority of functionally independent first integrals is installed. 

The same property of functional ambiguous of first integrals we have for systems of or- 
dinary differential equations [95, pp. 262 - 263], for total differential systems (see Subsec- 
tion 1.2), for linear homogeneous partial differential equations [52, p. 16], and for linear 
homogeneous systems of partial differential equations (Theorem 2.1). 

Example 4.1. Consider the Pfaff system of equations 

io^{x) = 0, LO2ix) = 0, (4.7) 
where the linear differential forms 

LOi{x) = Xi{l + X2) dxi +X2(1 — X2) dx2 + (^3 + X2X^) dx^ + (X4 + X2X3) dx^ for all x G M.^, 
U2{x) = Xi dxi — X2 dx2 + X4 dx^ + X3 dx^^ for all x G M^. 

27 



V.N.Gorbuzov 



Integral equivalence of multidimensional differential systems 



We have 



2a;i(x) + 2(1 - X2)oj2{x) = d{2x\ + {x^ + x^f) for all x G M^, 
2u^{x) - 2(1 + X2)oJ2{x) = d{2xl + (xg - x^f) for all x G M^. 



Therefore, by Definition 4.1, the functions 

F^:x^2x\ + (x3 + x^f for all x G M^, 



(4.8) 




is a first integral on the space of the Pfaff system of equations (4.7). Indeed, the function 
F^{x) = (Fi(x) - F2{x))/2 for all x eR"^. 



Definition 4.3. A set of functionally independent first integrals on the domain G' C G 
of system (Pf) with G G{G), j = 1, . . . ,m, is called a basis of first integrals {integral 
basis) on the domain G' of system (Pf) if for any first integral ^' on the domain G' 
of system (Pf), we have '^{x) = ^{F{x)) for all x G G', where $ is some function of 
class C^(EF), EF is the range of the vector function (2.3). The number k is said to be the 
dimension of basis of first integrals on the domain G' of system (Pf). 

Prom Definition 4.3 and Property 4.1, we get the following 

Property 4.2. The scalar functions (4.5) are a basis of first integrals on a domain O C G 
for the Pfaff system of equations (Pf) with m = n. 

Example 4.2. The Pfaff system of equations 



has the nonsingular matrix (4.2) on the set E, = {x: x^ ^ 0} (the determinant of this matrix 
is det w{x) = X2 for all a; G S). By Property 4.2, restrictions of the functions F^: x ^ x^ 

for all X G M'^, ^ = 1,2,3, are a basis of first integrals on any domain C S of the Pfaff 
system of equations (4.11). 

Under the condition X2 = the Pfaff system of equations (4.11) is the first-order ordinary 
differential equation dxi + 2 dx^ = 0. This differential equation has the general integral 
F: {xi, X3) Xi + 2x3 for all (x^^ , X3) G M?. 

4.3. Existence criterion of first integral 

The Pfaff system of equations (Pf) induces the linear differential forms (0.2). We add 
n — m linear differential forms 



iOi{x) = 0, L02{x) = 0, IJ^six) = 



(4.11) 



with the linear differential forms 



uJi{x) = dxi + dx2 + 2 dx^ for all x G M^, 

iL>2ix) = dx^ + 2dx2 + 2dx^ for all x G M^, 
u;3(x) = dxi + dx2 + (2 + X2) dx^ for all x G 



n 




for all X G G, = m + 1, . . . ,n. 



(4.12) 



i=l 
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with coefficients w^- G C(G), C = m + 1, . . . , n, i = I, . . . ,n, to the 1-forms (0.2) such that 
the set of the hnear differential forms (4.12) and (0.2) 

n 

uj^{x) = "w^iix) c?Xj for all X G G, ^ = 1, . . . ,n, (4-13) 

1=1 

are not linearly bound on the domain G. We form the square matrix of order n 

w{x) = \\w^i{x)\\ for all x G G. (4.14) 

The matrix (4.14) is nonsingular on a domain O C G with jiCc^ = 0. Then the matrix 
(4.14) on the domain Q, has the inverse matrix 

gr(a;) = ||g,.^(a;)|| for all X G fi. (4.15) 

The matrix (4.15) is a nonsingular on the domain square matrix of order n and 

w{x)g{x) = E for all x G (4.16) 

where E is the identity matrix of order n. 

We build the n linear differential operators of first order 

fi 

= 9d^)dx^ for all X G i = l,..., n, (4.17) 

which are not linearly bound on the domain J7 (because the matrix (4.15) is nonsingular on 
the domain Q). 

The operators (4.17) and 1-forms (4.13) are called contragredient if the coordinate relations 
(4.16) are hold. 

By the contragredient operators (4.17) and 1-forms (4.13), using the identity (4.16), we 
get the total differential of any scalar function F G C^{Q,) have the form 

n 

dF{x) = Y (5iF{x)i0i{x) forallxGi^. (4.18) 

i=l 

By virtue of (4.18) and Definition 4.1, we obtain an existence criterion of a first integral 
for a Pfaff system of equations. 

Theorem 4.2. A scalar function F G G^(J7) is a first integral on a domain Q C G of 
the Pfaff system of equations (Pf) with lVj G G(G), j = 1, . . . , m, if and only if the following 
conditions hold 

0^F(x) = O for all X GO, C = m + l,...,n. (4.19) 

4.4. Integral equivalence with linear homogeneous system of partial differential 
equations 

By the existence criterion of a first integral for the Pfaff system of equations (Pf ) (Theorem 
4.2), using a linear homogeneous system of partial differential equations, we can build an 
integral basis of the PfafF system of equations (Pf). 

Theorem 4.3. The scalar functions (2.2) are a basis of first integrals on a domain G'cG 
for the Pfaff system of equations (Pf) with ojj G C(G), j = 1, . . . ,m, if and only if the 
functions (2.2) are a basis of first integrals on the domain G' C O C G for the linear homo- 
geneous system of partial differential equations 

(S^{x)y = 0, C = m + l,...,n, (4.20) 

induced by the operators (4.17). 
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Proof. From the system of identities (4.19), we get the function F: Q — > M is a first integral 
on a domain ft C G oi the hnear homogeneous system of partial differential equations (4.20) 
(by Definition 2.1). Then, by the definitions of integral bases for the Pfaff system of equa- 
tions (Definition 4.3) and for the linear homogeneous system of partial differential equations 
(Definition 2.2), from Theorem 4.2, we obtain the criterion formulated in Theorem 4.3. ■ 

The systems (Pf) and (4.20) are called contragredient. 

Definition 4.4. We'll say that a Pfaff system of equations and a linear homogeneous 
system of partial differential equations are integrally equivalent on some domain if on this 
domain each first integral of the first system is a first integral of the second system and on the 
contrary each first integral of the second system is a first integral of the first system. 

By Theorem 4.3, the linear homogeneous system of partial differential equations (4.20) is 
integrally equivalent on the domain to the contragredient Pfaff system of equations (Pf). 

We supplement the l-forms (0.2) to the 1-forms (4.13) with only one condition to the 
linear differential forms (4.12): the 1-forms (4.13) is not linearly bound on the domain ft. At 
this viewpoint, the contragredient linear homogeneous system of partial differential equations 
(4.20) to the Pfaff system of equations (Pf) is constracted ambiguously. At the same time 
it does not influence (accurate within functional expression of basis integrals) on an integral 
basis of the Pfaff system of equations (Pf) and is regulated by Theorem 4.1. 

In regard to the domain Q, we have this domain is established by the domain of definition 
G of the Pfaff system of equations (Pf) and corrected by the possibility of construction of 
the inverse matrix (4.15) to the matrix (4.14). 

Example 4.3. Let us consider the Pfaff system of equations 

uj^{x) = 0, Lo^ix) = (4.21) 

with the 1-forms 

uJi{x) = dxi — dx2 — iXiX2 + X2 — 2,x| — 2x3X4) X2^{x^ — x^)~^ dxg -|- 
+ {X1X2 + X2 — 2x3X4 — 2x1)x2^{x3 — x^)'^ dx^ for all x G S, 

0^2(3;) = dxi + dx2 — {X1X2 — X2 + 2x| -I- 2x3X4) X2^ {x^ — X4)~^ dx^ + 

+ {X1X2 — X2 + 2x3X4 -|- 2x4)x^"^(x3 — X4)~^ dx4 for all x G S, 

where H = {x: X2 7^ 0, X4 7^ X3}. 
We add two 1-forms 

^al^) — ^3^2 ^(^3 ~ ^4)"^ "^^3 ~ 2:4X2 ^(3^3 — X4)~^ dx4 for all x G S, 

'^a{x) = — {x^ — X4)~^ dx^ + (X3 — X4)~^ dx4 for all x G S 

to the linear differential forms o;^ and uj2- 

The square matrix of fourth order (4.14) is generated by the coefficients of the 1-forms 
cjj , i = 1, . . . ,4. Since the determinant det'u;(x) = 2x^''^(x3 — X4)~^ 7^ for all x G S, we 
see that the matrix w is nonsingular on the set S. Therefore the 1-forms lo^, i = 1, . . . ,4, 
are not linearly bound on any domain $7 C S. 

We introduce the linear differential operators 

01 (x) = 0,5 d^^ - 0,5 d^^ for all x G S, 62 (^c) = 0,5 d^^ + 0,5 d^^ for all x G S, 

03(x) = 2(x3 + X4) d^^ + X2 d^^ + X2 d^^ for all x G S, 

(84 (x) = - Xi d^^ + X2 d^^ + X4 d^^ + X3 d^^ for all x G S, 

which are contragredient to the 1-forms co^, i = 1, . . . ,4. 
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The linear homogeneous system of partial differential equations 

03(a;)y = O, <3^{x)y = (4.22) 

is contragredient to the Pfaff system of equations (4.21). 
The restrictions of the functions 

F-^: X ^ Xi{x^ — x^)~^ for all x G {x: X4 7^ X3}, 

(4.23) 

F^-.x^xj {xl - (X3 + X4)2) for all x G M"^ 

are a basis of first integrals [96, p. 200; 34; 41] on the domain J7 of system (4.22). 

By Theorem 4.3, the restrictions of the functions (4.23) are an integral basis on the domain 
O of the Pfaff system of equations (4.21). 

4.5. Transformation of a Pfaff system of equations by known first integrals 

Theorem 4.4. // the Pfaff system of equations (Pf) has k, 1 ^ k ^ m, functionally 
independent on a domain G' <Z G first integrals (2.2), then this system on a domain Q. <Z G' 
with fxCo'^ = can be reduced to the form 

dF^{x) = 0, ^=l,...,k, 

(4.24) 

'^j^W = 0> Ja ^ {!'••• X = l,...,m-k. 

by a nonsingular linear transformation of the 1-forms (0.2). 

Proof. If the functions (2.2) are first integrals on a domain G' C G of system (Pf), then, 
by Definition 4.1, there exist the scalar functions b^^ G G^{G), ^ = 1, . . . ,k, j = 1, . . . ,m, 
such that the total differentials 

m 

dF^{x) = ^ b^j{x)uj{x) for all x G G', ^ = l,...,k. (4.25) 

Since the first integrals (2.2) are functionally independent on a domain G' and the 1-forms 
(0.2) are not linearly bound on the domain G, we see that the matrix 6(x) = ||^^j(a;)||^^ for 

all X £ G' (this matrix is induced by the coefHcients of expansion (4.25)) has rank6(x) = k 
for all X G rJ, where Q is a domain such that O C C and fiCc'^ = 0. 

Without loss of generality it can be assumed that the square matrix 6(x) = ||^^j(3^)|| for 

all X E Cl of order k (we obtain the matrix b from the restriction on the domain J7 of 
the matrix b by deletion of the last m — k columns) is nonsingular on the domain $7 (it 
always can be received by renumbering the 1-forms lUj , j = l,...,m). Then the system 
(Pf) is transformed into the system (4.24) under the nonsingular on the domain linear 
transformation of the 1-forms 

rn 

l^ = ^b^j{x)ujj, ^ = l,...,k, l0=uj0, e = k+l,...,m. ■ 

The differential system (4.24) is constructed with the help of the not linearly bound on 
the domain G differential forms 

n 

'^^i^) — 9^ F^{x) dxj for all X G ^ = 1, . . . , k, 

i=l 

^j^ix) for all X G ri, j^^ G {1, . . . ,m}, X = 1, . . . ,m — k. 

In this connection, we have 
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Theorem 4.5. The restrictions on a domain 0, c G' , /iCg'^ = 0, of the functions (2.2) 
are k, 1 ^ k ^ m, functionally independent first integrals on the domain J7 of the Pfaff 
system of equations (4.24). 

By Theorems 4.4 and 4.5, we obtain an existence criterion of functionally independent 
first integrals for a Pfaff system of equations. 

Theorem 4.6. The Pfaff system of equations (Pf) has k, 1 ^ k ^ m, functionally inde- 
pendent first integrals on a domain Q G G if and only if this system on the domain J7 can 
be reduced to the form (4.24) by the nonsingular linear transformation of the 1-forms (0.2). 

4.6. Closed systems 

Definition 4.5. The Pfaff system of equations (Pf) is called closed on a domain (Z G if 
a contragredient linear homogeneous system of partial differential equations (4.20) is complete 
on the domain Q.. 

For example, the contragredient linear homogeneous system of partial differential equati- 
ons (4.22) to the Pfaff system of equations (4.21) is complete on a domain C S. Indeed, 
since the basis of first integrals (4.23) has the dimension n — m = 4 — 2 = 2, we see that 
the contragredient system (4.22) is complete (by Theorem 2.13). Thus the Pfaff system of 
equations (4.21) is closed on the domain $7. 

Theorem 4.7. The Pfaff system of equations (Pf) is closed on a domain H <Z G if and 
only if this system on the domain H has an integral basis of dimension m. 

Proof. By Theorem 2.10, a basis of first integrals of the complete system (4.20) on a 
neighbourhood of any point of its normalization domain H has the dimension m. Taking 
into account Theorem 4.4, we obtain the closed Pfaff system of equations (Pf) and the con- 
tragredient complete linear homogeneous system of partial differential equations (4.20) have 
the same dimensions of integral bases. Therefore these dimensions are equal m. 

Thus, by Theorem 2.13, the system (Pf) is closed on a domain H d G \i and only if an 
integral basis of this system on the domain H has the dimension m. ■ 

Theorem 4.7 is a criterion of closure for a Pfaff system of equations in the terms of the 
dimension of an integral basis. 

Example 4.4. The Pfaff system of equations (4.7) has two equations (m = 2) and two 
functionally independent on the space first integrals (4.8) and (4.9). 

Therefore the system (4.7) is closed and the functions (4.8) and (4.9) are an integral basis 
of system (4.7) on the space R^. 

Likewise, by Theorem 4.7, we can prove that the Pfaff system of equations (4.11) (see 
Example 4.2) is closed on space M^. Indeed, this system has an integral basis of the dimension 
three on space M^. 

In Theorem 4.7 we can take H as a normalisation domain (Definition 2.4) of the contra- 
gredient system (4.20) to the system (Pf). 

Prom Theorem 4.6 under the condition k = m and Theorem 4.7, we obtain the following 
criterion of closure for a Pfaff system of equations [53, pp. 110 - 111]. 

Theorem 4.8. The Pfaff system of equations (Pf ) is closed on a domain <Z G if and 
only if this system on the domain J7 can be reduced to the differential system 



by the nonsingular linear transformation of the 1-forms (0.2). 

The scalar functions F- G C^(0), j = l,...,m, are first integrals on the domain Q, 
both for the system (4.26) and the system (Pf). The systems (Pf) and (4.26) are integrally 
equivalent on the domain Q. d G. 

Example 4.5. The Pfaff system of equations (4.7) is transformed into the system 



dFAx) = Q, i = l,... 



m 



(4.26) 



??i(x)=0, 772(0;) = 0, 



(4.27) 



where the 1-forms 
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r]i(x) = 4:X^ dxi + 2(^3 + x^) dx^ + 2{x^ + x^) dx^ for all x G M^, 

%(^) — ^2 + 2(^3 ~ ^4) ^^^3 ~ 2(a;3 — a;4) da;4 for all x G R^, 
under the nonsingular on the space linear transformation of the 1-forms u-^ and 0J2 
r/i = + 2(1 - X2)uj2 , r?2 = - 2(1 + X2)w2. 

Since 

r?i(x) = c?(2xf + (X3 + X4)2) for all x G M^, 

7/2(2;) = d{2xl + (0:3 - X4)2) for all x G M^, 
we see that the system (4.27) can be reduced to the form 

d{2xl + (X3 + x^f) = 0, d{2xl + (X3 - x^f) = 0. 

By Theorem 4.8, the system (4.7) is closed and the first integrals (4.8) and (4.9) are 
an integral basis on the space of system (4.7) (using Theorem 4.7, the same result was 
obtained in Example 4.4). 

Using the nonsingular on the space linear transformation of the 1-forms oj^ and (jJ2 

CTi = 2UJ^ + 2(1 - X2)UJ2 ' <^2=^2^ 

we get the system (4.7) can be reduced to the system 

(Ti(x) = 0, (T2(a;) = 0, (4.28) 

where the 1-forms 

(T-i{x) = Axi dxi + 2(^3 + x^) dx^ + 2{x^ + X4) dx^^ for all x G M^, 

o"2(x) = a:;;^ dx^ — X2 dx2 + ^4 dx:^ + X3 dx^^ for all x G M^. 

Since 

0-1 (x) = d{2x\ + (X3 + X4)2) for all x G M^, 
2(72 (x) = d{x\ — X2 + 2x3X4) for all x G M^, 
we see that the system (4.28) can be reduced to the form 

d{2xl + (X3 + X4)2) = 0, d{xl -xl + 2x3X4) = 0. 

By Theorem 4.8, the system (4.7) is closed and the first integrals (4.8) and (4.10) are an 
integral basis on the space of system (4.7). 

4.7. Interpretation of closure in terms of differential forms 

In [58], the interpretation of complete solvability for the total differential system (TD) 
in terms of differential forms was given. We give the interpretation of closure for the Pfaff 
system of equations (Pf) in terms of differential forms. 

Lemma 4.1. Suppose the linear differential forms ojp G C°°{G), p = 1, . . . , s. Then the 
system of exterior differential identities 

dwp{x) ^OJ'^{x) ^ . . . /\u)g{x) = for all X e G, p = l,...,s, (4.29) 

is invariant under the nonsingular on the domain G linear transformation of the 1-forms 
ujp, p = l,.. .,s. 

Proof. Let the 1-forms , 6 = 1, . . . , s, be nonsingular on the domain G linear combi- 
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nations of the 1-forms , p = I, . . . , s, i.e., 

s 

hi^) = '^Spix)^pix) for all X G G, S=l,...,s, (4.30) 
p=i 

where the scalar functions "^g^ G C°°{G), d = l,...,s, p = l,...,s, the square matrix 
^{x) = for all X E G of order s is nonsingular on the domain G. Then the 

exterior product 

A IJx) = det^'(a;)f A ujJxj) for all x eG. 
p=i ^ V p=i ^ / 

Therefore the system of identities 

dlg{x) Al^{x) A ... Alg{x) = forallxGG, 6=1,. ..,s, (4.31) 
is valid if and only if 

dlg{x) Au}^{x) A . . . Au}g{x) = for all x G G, 5 = 1,. ..,s. (4.32) 
By the representations (4.30), the exterior differentials 

s s 

dk{x) = ^5p{x) dujp{x) + d^spi^) ^ ^p{x) for all X € G, 6 = 1,..., s. 
p=i p=i 

Thus the exterior products 

s 

dlgix) A Ui{x) A ... A a;,(x) = '^Sp{x)dUp{x) A u^ix) A ... A u,{x) 

P=i 

(4.33) 

for all X G G, 6=1,..., s. 

Using the identities (4.29), from the identities (4.33), we get the identities (4.32). Therefore 
the identities (4.31) are consistent. This yields that the system of identities (4.29) is invariant 
under the nonsingular on the domain G transformation (4.30) of the 1-forms Up, p= 1, . . . , s.M 

Theorem 4.9. If the Pfaff system of equations (Pf) with cOj G G°°(G), j = 1, . . . ,m, has 
m functionally independent on a domain G' C G first integrals, then the exterior products 

dwj{x) Auj-i^{x) A . . . Aujjnix) = for all x G G', j = l,...,m. (4.34) 

Proof. Let the Pfaff system of equations (Pf) has the m functionally independent on a 
domain G' C G first integrals 

Ff.x^ F-{x) for all x G G', j = l,...,m. (4.35) 

By Definition 4.1, the total differentials 

m 

dFj{x) = Y bjQ{x)uj^{x) for all x G G', j = l,...,m, (4.36) 
C=i 

where the scalar functions b^^ G C°°{G'), j = 1, ... ,171, = 1, . . . ,m, the square matrix 
6(x) = ||&j^(x)|| for all x G G' of order m is nonsingular on a domain CI C G' with 
pCqiQ = (see the proof of Theorem 4.4). 

By the Poincare theorem [97, p. Ill] (for any differential g-form 6 G G°°(G), we have 
the identity d{dQ{x)) = for all x G G), in view of the first integrals (4.35) of system (Pf) 
we obtain 
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(fF-{x) = for all x G G', j = 1, . . . , m. (4.37) 

For the 1-forms 

m 

lj{x) bji^{x)uj^{x) for all x e G', j = l,...,m, (4.38) 

C=i 

using the identities (4.36) and (4.37), we get the exterior products 

dlj{x) Mj{x) /\...A l^{x) = for all x G G', j = l,...,m. (4.39) 

Since the matrix b is nonsingular on a domain C G', we see that from the represen- 
tations (4.38) it follows that the 1-forms 1^ , j = 1, . . . ,m, on the domain J7 are the result 
of the nonsingular on the domain J7 linear transformation of the 1-forms ujj , j = 1, . . . ,m. 
Therefore, by Lemma 4.1, the identities (4.39) are valid if and only if the identities hold 

dLOj{x) A ujj{x) A ... A uJ^{x) =0 for all x G j = 1, . . . , m. 

Whence, using uoj € C°°{G), j = 1, . . . ,m, fiCcft = 0, we get the identities (4.34). ■ 
Prom Theorems 4.7 and 4.9 we obtain the following 

Theorem 4.10. // the Pfaff system of equations (Pf) with oj^ G G°°(G), j = 1, . . . ,m, 
is closed on a domain H G G, then the exterior products 

dLOj{x) A iOi{x) A . . . A iOjj^(x) =0 for all x G if, j = 1, . . . , m. 

Theorem 4.11. // the system of the exterior identities 

duj{x) AuJi{x) A . . . Au^{x) = for all x e G, j = l,...,m, (4.40) 

is valid, then the Pfaff system of equations (Pf) with coj G G°°(G), j = 1, . . . ,m, is closed 
on a domain O C G, where Q is a domain such that fi Cq^ = 0. 

Proof. Under the condition n — m = 0, we have the exterior products 

dLOj{x) A iOi{x) A . . . A C0jjj_{x) for all x G G, j = 1, . . . ,m, 

are (n -|- 2) -forms of n variables and under the condition n — m = 1, we get these exterior 
products are (n -|- l)-forms of n variables. In these cases the identities (4.40) are valid on 
the domain G. Now let us prove that the system (Pf ) with n — m = and the system (Pf ) 
with n — m = 1 are closed on a domain $7 C G, CgO = 0. 

If n — m = 0, then the Pfaff system of equations (Pf) on the domain 0, has n first 
integrals (4.5) (by Property 4.1). These first integrals are an integral basis of system (Pf) 
on the domain (by Property 4.2). Taking into account Theorem 4.7, we obtain the Pfaff 
system of equations (Pf) with n — m = is closed on the domain J7. 

Suppose n — m = 1. Then the matrix (4.2) of system (Pf) is an (n — 1) xn matrix and has 
Tankw{x) = n—1 for all x G ^l, where O C G is a domain such that jJ-Cc^ = 0. Therefore 
the system (Pf ) on some domain 17 C $7 with /i C^^ft = can be solved for m = n ~ 1 
differentials. For example, if tlie system (Pf) is solved for dxi,. . . ,dx^_i , then the system 
(Pf) can be reduced to the system of n — 1 ordinary differential equations 

dx 

^ = P^(xi,...,xJ, T = l,...,n-1, (4.41) 

with the right hand sides G G°°(0), r = 1, . . . , n — 1. 

The system (4.41) has n — 1 functionally independent on the domain ft first integrals (by 

Theorem 1.3 with m = 1). Hence the system (Pf) with n — m = 1 has an integral basis of 
dimension m = n — 1 on the domain fl, where O is a domain from G such that /i Cq^I = 
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(because IJ'C^i} = and /iCg^ = 0). By Theorem 4.7, the Pfaff system of equations (Pf) 
with n — m = 1 is closed on the domain $7. 

Thus the assertion of Theorem 4.11 is vahd for the Pfaff system of equations (Pf) of 
codimension null (n — m = 0) and codimension one (n — m = l). 

The proof of Theorem 4.11 for the Pfaff system of equations (Pf) of codimension s > 1 
{n — m = s) is by mathematical induction on s. 

We assume that Theorem 4.11 is true for s > 1, s = n — m, i.e., if the system of exterior 
differential identities (4.40) with m, = n — s, s > 1, is valid, then the Pfaff system of equations 
(Pf ) of codimension s is closed on a domain O C G, Cq^ = 0. 

By Theorem 4.8, we assume that if the system of exterior differential identities (4.40) 
with m = n — s, ,s > 0, is valid, then there exists a nonsingular on the domain Q C G, 
hCg^ = 0, linear transformation of the 1-forms ojj , j = l,...,m, such that the Pfaff 
system of equations (Pf) of codimension s can be reduced to the differential system 

dF^{x) = 0, j = 1, . . . ,m, m = n — s, s > 1. 

Let us consider the Pfaff system of equations (Pf) of codimension s + 1. 

Let the system of exterior differential identities (4.40) with m = n— (s+1), s > 1, be valid. 
If we fix , then the identities (4.40) with m = n — (s + 1), s > 1, are corresponding to the 
identities (4.40) with m = n — s, s > 1. Therefore, by the inductive assumption, there exists 
a linear transformation of the 1-forms ujj , j = 1, . . . ,m, on the domain $7 C G, juCcO = 0, 
such that the Pfaff system of equations (Pf) with m = n— (s + 1), s > 1, can be reduced to 
the Pfaff system of equations 

lj{x)=0, j = l,...,m, m = n-{s + l), s>l, (4.42) 

with the 1-forms 

lj{x) = dFj{x) + Qjix) dx„ for all x G fi, j = 1, . . . , n — s — 1, s > 1. (4.43) 

Using the Poincare identities (4.37) with m = n — (s -|- 1), s > 1, we obtain the exterior 
differentials of the 1-forms (4.43) 

dlj{x) = dgj{x) A dx^ for all x e 0,, j = 1, . . . ,n — s — 1, s > 1. (4.44) 

On the other hand, the 1-forms Ij , j = l,...,n — s — 1, s > 1, are the result of 
the nonsingular on the domain Q C G, n CcO = 0, linear transformation of the 1-forms 
ujj, j = 1, . . . ,n — s — 1, s > 1. Then, by Lemma 4.1, the exterior products 

dlj{x) A li{x) A . . . A Ijnix) =0 for all x G j = 1, . . . ,n — s — 1, s > 1. 

Therefore the exterior differentials 

m 

dlj{x) = ^ Qj(^{x) A /^(x) for all X G Jl, j = 1, . . . ,n — s — 1, s > 1. (4.45) 
C=i 

Combining the identities (4.44) and (4.45), we obtain the total differentials 

m m 

d9j{x) = hj^{x){dF^{x) + g^{x) dx^) = ^ /i^-^(x) dF^{x) + hj{x) dx^ 
C=i C=i 

for all x G Q, j = 1, . . . ,n — s — 1, s > 1, 

where 

m 

hj{x) = hji^{x)g^{x) for all X G il, _7 = 1, . . . , n — s — 1, s > 1. 

C=i 
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From these identities it follows that 

= 9j{Fi{x), . . .,F^{x),Xn) for all x e n, j = 1, . . . ,n - s - 1, s > 1, 

where g^, j = 1, . . . ,n — s — 1, s>l, are holomorphic scalar functions of m + 1 variables 
-^1 ' ■ ■ ■ ' jdx^ . 

Then the equations (4.42) with the 1-forms (4.43) are 

dPj + 9j{Fi,---,F^,x^)dx^ = ^, j = 1, s > 1. 

This system is a system of m ordinary differential equations and has a basis of first 
integrals of dimension m, m = n — (s + l), s > 1 (by Theorem 1.3 with m = 1). Moreover, 
this basis is an integral basis on the domain Q of system (Pf) with n — m = s + 1, s > 1. 

By Theorem 4.7, the system (Pf) with n — m = s + 1, s > 1, is closed on the domain il. ■ 

From the proof of Theorem 4.11, we get the following statements. 

Corollary 4.1. The Pfaff system of equations (Pf) with Uj G C°°{G), j = 1, . . . ,m, and 
m = n — 1 is closed on a domain C G with I-lCg^ = 0. 

Corollary 4.2. The Pfaff system of equations (Pf) with lJj G C°°(G), j = 1, ...,m, 
and m = n — 1 can be reduced to the integrally equivalent on a domain Q G G, fx CgJI = 0, 
system of n ~ 1 ordinary differential equations (4.41) by the nonsingular on the domain O 
linear transformation of the 1-forms (0.2). 

Using Theorems 4.10 and 4.11, we obtain the Frobenius theorem [53, pp. 110 - 112; 97, 
pp. 131 - 136], which is a criterion of closure for a Pfaff system of equations with the help of 
exterior products of differential forms. 

Theorem 4.12. The Pfaff system of equations (Pf) with coj G C°°{G), j = l,...,m, 
is closed on a domain G G, fiCo^ = 0, if and only if the system of exterior differential 
identities (4.40) is valid. 

The system of exterior differential identities (4.40) is called [58, p. 302] the Frobenius 
conditions of closure for the Pfaff system of equations (Pf). 

Example 4.6. The Pfaffian differential equation 

yz dx + 2xz dy + 3xy dz = (4.46) 

induces the vector field A: {x,y,z) — > {yz,2xz,3xy) for all {x,y,z) G M'^ with the rotor 
rot A: {x,y, z) ^ {x, — 2y, z) for all (x, y, z) G M^. 

The scalar product A{x,y, z) rot A{x,y, z) = for all {x,y, z) gM.^, i.e., the vector field A 
is orthogonal to the rotor of A. This condition is equivalent to the Frobenius condition (4.40). 

Therefore the Pfaffian differential equation (4.46) is closed on space and this equation 
has an integral basis of dimension one. 

The first integral 

F: {x,y,z) — >■ xy'^z^ for all {x,y,z) G 
is an integral basis of the Pfaffian differential equation (4.46). 
4.8. Nonclosed systems 

Let us consider the Pfaff system of equations (Pf ) such that the contragredient linear ho- 
mogeneous system of partial differential equations (4.20) is incomplete on a domain O C G. 
In this case, the Pfaff system of equations (Pf) is said to be nonclosed on the domain O. 

Further, adding the equations of the forms (2.9) to the system (4.20), we get a correspon- 
ding complete system to the incomplete system (4.20) and the defect 6, 1 < 6 ^ m, of system 
(4.20). For this complete system we obtain a normalization domain H G ^ (Definition 2.4). 
Then, by Theorem 4.3, we have 

37 



V.N.Gorbuzov 



Integral equivalence of multidimensional differential systems 



Theorem 4.13. The nondosed on a domain C G Pfaff system of equations (Pf) 
has a basis of first integrals of dimension m — S on a normalization domain H G ^ of the 
contragredient linear homogeneous system of partial differential equations (4.20), where S is 
the defect of system (4.20). 

Let us remember that the complete system (4.20) has the defect 6 = 0. Then, using 
Theorems 4.7 and 4.13, we obtain the generahzing statement about a basis of first integrals 
for the closed or nonclosed Pfaff system of equations (Pf). 

Theorem 4.14. The Pfaff system of equations (Pf) on a normalization domain H C G 
of the contragredient linear homogeneous system of partial differential equations (4.20) has a 
basis of first integrals of dimension m^6, where 5 is the defect of system (4.20), ^ 5 ^ m. 

Example 4.7. Consider the Pfaff system of equations 

Ui{x) = 0, uj2{x) = (4.47) 

with the 1-forms 

uji{x) = dxi + dx2 + dx'^ + dXi^ for all a; G M^, 

u)2{x) = dxi + 2dx2 + x^dx^ + dx^ for all x G M^. 

We add two 1-forms 

0J^{x) = dx^ for all x G R^, 0^4 (x) = dx^^ for all x G M^. 

to the linear differential forms Ui and 0^2 ■ 

The linear differential forms u^, i = 1, . . . ,4, are not linearly bound on the space R^. 
Using the not linearly bound on the space contragredient linear differential operators 

©i(x) = 2d^^ - d^^ for all x G R^, 62 (x) = -d^^+ d^^ for all x G R^, 

03(x) = (X4 - 2)d,^ + (1 - x^)d,^ + d,^ for all x G R^ 

04 (x) = -d^^+ d^^ for all x G R^ 

to the 1-forms , i = 1, . . . ,4, we obtain the contragredient linear homogeneous system of 
partial differential equations 

03(x)y = O, 04(x)y = O (4.48) 

to the Pfaff system of equations (4.47). 
Since the Poisson bracket 

= [<S>^{x), = 9., = <»2(^) for all x G R^ 

we see that the system (4.48) is incomplete. Therefore the Pfaff system of equations (4.47) is 
nonclosed. 

The system (4.48) with the help of the operator ©34 can be reduced to the complete 
system 

03(x) y = 0, 04(x) y = 0, ©34 (x) y = 0, 

The first integral 

F: X ^ X;^ + ^2 + X3 + X4 for all x G R^ 

is an integral basis on the space R^ of this complete system. 

This function is a basis of first integrals on the space R^ of the nonclosed Pfaff system of 
equations (4.47). 
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4.9. Integral equivalence with total differential system 

One more approach for building of an integral basis of a PfafF system of equations is based 
on a reducing this system to integrally equivalent total differential system. 

Definition 4.6. We 'U say that a Pfaff system of equations and a system of total differential 
equations are integrally equivalent on some domain if on this domain each first integral of 
the first system is a first integral of the second system and on the contrary each first integral 
of the second system is a first integral of the first system. 

The linear differential forms (0.2) are not linearly bound on the domain G. Therefore the 
mxn matrix (4.2) has rank'u;(a;) = m for all x G where a domain O C G and /xCqO = 0. 
Then, the square matrix w{x) = \\wj^{x)\\ for all x E G of order m is nonsingular on the 

domain Q. Thus the Pfaff system of equations (Pf) can be reduced to the system of total 
differential equations 

n 

dxj= ^ aj^{x)dx^, j = l,...,m. (4.49) 

i'=m+l 

The Pfaff system of equations (Pf) and the system of total differential equations (4.49) 
are integrally equivalent on some domain G' C O C G, i.e., we have the following assertions. 

Theorem 4.15. A scalar function F G G^(G') is a first integral on a domain G' C G of 
the Pfaff system of equations (Pf) with utj € C{G), j = 1, . . . ,m, if and only if this function 
is a fi,rst integral on the domain G' of the system of total differential equations (4.49). 

Theorem 4.16. The scalar functions (2.2) are a basis of first integrals on a domain 
G' G G of the Pfaff system of equations (Pf) with Uj € C{G), j = 1, . . . ,m, if and only if 
these functions are a basis of first integrals on the domain G' of the system of total differential 
equations (4.49). 

Using Theorems 3.6, 4.7, and 4.16, we can prove the following 

Theorem 4.17. The Pfaff system of equations (Pf) is closed on a domain G' C G if and 
only if the total differential system (4.49) on the domain G' is completely solvable. 

Example 4.8. The Pfaff system of equations 



(4.50) 



(4.51) 



2xi{l + X2) dxi + 6x2 dx2 + 3x3(2 + X2) dx^ + 3x4(2 + X4) dx^ = 0, 

4x;^(l + x^) dxi — 6x2 dx2 + 3x3(1 + 2xi) dx^ + 3x4(1 + 2xj^) dx^ = 

can be reduced to the system of total differential equations 

_ 3 _i 3 _i 

dx-^ — ~ X^X-^ dx^ ~ X^X-^ dx^ ^ 

dxi2 — ~ '^3'^ 2 dx^ ~ x^x<2 dx ^ J 

which is defined on the set H = {x: x^ 7^ 0, X2 ^ 0, 3 + 2xi + X2 7^ 0}. 
Since the Poisson bracket 

^X-^ " 2 ^3^1 ^ ^X^ ~ 2 ^3^2 ^ ' ^X2 ~ 2 ^^^^ ^ ^^3 ~ 2 ^^^^ ^ ^^4 ~ ^ 

for all X G S, S = {x: x^ 7^ 0, X2 7^ 0}, 

we see that the system (4.51) is completely solvable on any domain G C S. 

By Theorem 4.17, the Pfaff system of equations (4.50) is closed on any domain G C S. 
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The functionally independent first integrals 

Fi:x^2xl + 3x1 + ^^4 for all x G M^, 

F2: X ^2xl + xl + xl for all x G 

are a basis of first integrals on any domain G C H of the system (4.51), and a basis of first 
integrals on any domain G' C S of the system (4.50). 
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